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ABSTRACT 


This  thesis  is  a  continuation  of  previous  work  to  apply 
modern  multiparameter  estimation  techniques  to  the  problem 
of  estimating  attrition  rates  for  a  large  number  of  small 
inventory  cells  in  manpower  planning  models  used  by  the  U.S. 
Marine  Corps.  The  main  advances  involve  the  promising 
introduction  of  empirical  Bayes  (non-constant  shrinkage) 
techniques,  recognition  of  the  non  symmetric  nature  of  the 


errors  with  a  response  to  this,  and  some  insight  into  all 
aggregation  plans  that  should  help  provide  greater  stability 


for  the  estimation  methods.  In  addition,  the  roles  of  some 
middle  level  methodological  choices  are  explored.  (  \  -  ^  ^ 
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within  the  time  available,  to  ensure  that  the  programs  are 
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I.  INTRODUCTION 

A.  PURPOSE 

In  support  of  project  OPUS  (Officer  Planning  and  Utility 
System),  Headquarters,  U.S.  Marine  Corps,  in  1985,  requested 
assistance  in  exploring  new  methods  of  generating  manpower 
loss  rates  in  order  to  improve  upon  the  currently  used 
method  of  force  rate  analysis.  In  response,  students  and 
faculty  in  the  Operations  Research  department  at  the  Naval 
Postgraduate  School  began  exploring  the  use  of  modern 
multiparameter  estimation  techniques  for  this  problem. 
Special  emphasis  is  placed  on  the  small  cell  problem,  i.e., 
categories  of  officer  skill,  grade  and  length  of  service 
which  have  low  inventory  figures.  Historically,  rate 
estimators  for  small  cells  are  unstable  and  a  large  number 
of  these  cells  exist. 

This  paper  builds  upon  previous  work  on  this  problem. 
The  main  advances  are  to  append  new  measures  of 
effectiveness  as  requested  by  the  Navy  Personnel  Research 
and  Development  Center  (NPRDC) ,  introduce  a  class  of 
empirical  Bayes  estimators,  and  explore  the  effects  of  some 
middle  level  choices  in  applying  the  new  and  existing 
methods.  The  ultimate  goal  is  to  refine  the  techniques 
presented  here  in  order  to  validate  a  clear  policy  for 
predicting  loss  rates. 


B .  BACKGROUND 


For  an  introduction  into  Marine  Corps  policy  concerning 
manpower  planning,  the  reader  is  referred  to  a  thesis 
submitted  in  September  1985,  at  the  Naval  Postgraduate 
School,  by  Major  D.D.  Tucker  [Ref.  1],  who  presents  a 
detailed  background  into  the  Marine  Corps  officer  structure 
and  the  manpower  planning  process.  Tucker  also  provides 
basic  attrition  rate  theory  and  calculates  overall  attrition 
rates  in  several  different  formats  as  they  pertain  to  the 
Marine  Corps.  As  the  aggregation  rate  begins  to  grow 
smaller  with  refinement,  Tucker  illustrates  the  irregular 
behavior  of  losses  due  to  voluntary  attrition  and  the  small 
numbers  of  losses  obtained  in  such  aggregates.  Tucker  [Ref. 
l:p.  50]  introduces  the  James-Stein  technique  of  loss  rate 
estimation  as  a  shrinking  of  individual  cell  averages  toward 
a  grand  average  in  order  to  reduce  the  risk,  or  squared 
differences  of  forecast  and  actual  values,  the  goal  being  an 
improvement  over  the  classical  maximum  likelihood  estimator, 
referred  to  from  now  on  as  the  MLE.  Tucker  chose  the  ranks 
of  First  Lieutenant  and  Lieutenant  Colonel  for  evaluation 
using  military  occupational  ppecialty  (MOS)  groups  of  combat 
support  and  ground  comb*-  -  (three  MOS ' s  were  selected  for 
each  group) .  These  aggregation  schemes  were  carried  through 
other  studies  and  are  included  for  historical  purposes  in 
this  paper. 
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Another  thesis  submitted  at  the  Naval  Postgraduate 
School  in  March  1986,  by  Major  J.R.  Robinson  [Ref.  2] 
focused  on  a  technique  called  limited  translation  of  the 
James-Stein  estimator.  This  particular  model  attempts  to 
minimize  the  risk  mentioned  above  by  reducing  the  shrinkage 
of  rate  estimates  for  cells  which  fall  outside  a  certain 
range  of  values  centered  on  the  grand  mean.  It  was  thought 
that  shrinking  all  cells  by  the  same  rate  toward  this  grand 
mean  may  be  unwise  since  the  attrition  rates  for  those 
particular  cells,  those  farthest  from  the  grand  mean,  may 
contain  attributes  which  are  different  than  the  majority  of 
the  cells  in  the  aggregate.  Additionally,  Robinson’s  work 
carries  forward  the  study  done  by  Tucker  concerning  MLE  and 
James-Stein  estimators. 

These  feasibility  studies  were  limited  in  scope  due  to 
the  format  of  the  summary  date  tape  available  (acquired  from 
NPRDC) ,  which  allows  only  coarse  cell  definition.  A  cell  is 
defined  to  be  a  cross-classification  of  the  forty  military 
occupational  fields  (OF) ,  thirty-one  lengths  of  service 
(LOS),  and  ten  grades  for  a  total  of  12,400  categories  for 
manpower  planning  purposes  [Ref.  2:p.  10].  The  data  does 
not  distinguish  between  limited  duty  officers,  those 
officers  specifically  designated  for  limited  duty  within 
certain  MOS's,  and  unrestricted  officers  (regular  and 
reserve)  [Ref.  l:p.  20].  These  factors  dictated  a  broad 
aggregation  scheme  that  created  stable  cell  inventories  and 
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allowed  the  common  variance  assumption  crucial  to  their 
work.  But  some  of  the  application  models  require  more 
refined  information  along  those  lines  and  experiments  with 
finer  aggregation  levels,  necessary  in  some  real  cases, 
caused  this  assumption  to  be  violated. 

Tucker  and  Robinson  did  conclude  that  the  current  scheme 
of  attrition  rate  prediction  could  be  improved  upon  by  their 
methods,  although  no  dominant  scheme  was  uncovered.  Tucker 
and  Robinson  did  identify  problem  areas  in  their  studies 
that  required  additional  work.  Tucker  [Ref.  l:p.  71]  cited 
a  need  for  a  better  aggregation  method  to  produce  a  more 
homogeneous  attrition  behavior  in  all  cells.  Also,  small 
probabilities  of  loss  within  a  cell  were  not  dealt  with 
successfully.  Robinson  [Ref.  2:p.  32]  revealed  the  small 
cell  problem  and  the  inability  to  normalize  the  cell  means 
or  stabilize  the  variance  with  a  data  transformation.  The 
recommendations  of  the  previous  studies  highlighted  the  need 
for  a  more  refined  data  tape  with  current  information  to 
include,  among  others,  full  Military  Occupational  Specialty 
(MOS)  information,  grade  separation  to  include 
regular/reserve  status,  promotion  zone  data  and  breakdown  of 
attritions  by  type. 

C.  PROGRESS 

The  project  of  loss  rate  estimation,  sponsored  by  NPRDC , 
is  currently  moving  into  the  implementation  stage.  An 
operational  data  tape  is  now  available  which  includes 
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detailed  information  on  Marine  Corps  officers  for  the  years 
1977-86.  Work  has  begun  at  the  Naval  Postgraduate  School  to 
break  out  the  data  in  a  useful  format  and  its  availability 
in  now  nearing  completion. 

In  addition  to  the  results  stated  earlier.  Tucker  and 
Robinson  provided  loss  rate  estimates  for  small  cells  with 
no  attrition,  i.e.,  MLE  equal  to  zero,  and  risk  effects  of 
the  different  schemes  studied.  These  results  provided  the 
springboard  for  this  study  and  the  requirements  for 
additional  future  studies. 

Captain  R.W.  Larsen  [Ref.  3]  has  developed  a  promising 
aggregation  scheme  based  on  cluster  analysis,  i.e.,  a 
classification  scheme  to  aggregate  cells  which  reflect  a 
greater  degree  of  homogeneity  in  attrition  rates  than  are 

i 

allowed  using  current  aggregation  methods.  Table  1  displays 
the  resulting  scheme,  given  in  year  of  current  service  (YCS) 
groups  proposed  for  implementation  by  Larsen. 


Table  1.  AGGREGATION  METHOD  PROPOSED  BY  LARSEN 


MOS  Category 

Bounded  YCS  Groups 

Fixed-Wing  Pilots 

(1-6,8-19) 

(7) 

(20-25) 

(26) 

Rotary- Wing  Pilots 

(1-5,8-19) 

(6,7) 

(20-25) 

(26) 

Naval  Flight  OfTicers 

(1-5,8-19) 

(6,7) 

(20-25) 

(26) 

Lawyers 

(1-6,8-19) 

(7) 

(20-25) 

(26) 

All  Else  • 

(1-3,6-19) 

(4,5) 

(20-25) 

(26) 
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The  parentheses  encompass  YCS  groups  which  behave  similarly. 
In  each  case  the  grade  is  fixed,  i.e.,  one  should  aggregate 
over  YCS  before  aggregating  over  grade. 

This  information  will  allow  for  more  successful 
application  of  the  loss  rate  estimation  techniques  of 
interest  here  when  the  new  data  tape  is  provided. 


D .  GOALS 

In  an  attempt  to  refine  the  methods  of  attrition  rate 
estimation  presente  1  in  the  previous  pilot  studies,  the 
following  goals  were  set  for  this  paper: 

1.  Modify  the  existing  methods  to: 

a)  use  the  refined  data  format, 

b)  extend  the  work  to  the  case  of  unequal  or  non 
constant  variance  in  the  small  cell  inventory 
problem, 

c)  study  alternate  transformation  inversion 
techniques , 

d)  introduce  and  evaluate  additional  measures  of 
effectiveness  suggested  by  NPRDC  to  include,  but 
not  limited  to,  cross  validation  underages  and 
overages  of  the  deviation  between  forecasts  and 
actuals , 

e)  examine  the  graphical  effect  of  shrinkage  in  the 
original  scale. 

2.  Introduce  empirical  Bayes  estimation  method  in  several 
forms  for  attrition  rates  for  consideration  as  another 
option  to  solve  the  problem. 
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II.  EXTENSIONS  TO  ESTIMATION  METHODS 


A .  GENERAL 

As  stated  in  Chapter  I,  Tucker  and  Robinson  published 
pilot  studies  introducing  the  James-Stein,  limited 
translation  James-Stein,  Maximum  Likelihood  and  Transform 
Scale  Cell  Average  estimators  as  alternatives  to  the 
aggregate  methods  previously  used  for  manpower  planning  in 
the  Marine  Corps.  Their  performance  was  promising,  and  this 
chapter  will  investigate  ways  to  enhance  their  performance 
and  possibly  develop  a  dominant  estimator  for  future 
attrition  prediction.  Additionally,  NPRDC  has  expressed 
interest  in  additional  measures  of  effectiveness  for  these 
estimators  which  may  lead  to  a  sharper  direction  for  the 
goal  of  producing  a  viable  estimation  policy.  This  chapter 
concentrates  on  the  following  issues  concerning  previous 
studies . 

B.  NON  CONSTANT  VARIANCE 

For  a  cell  having  an  inventory  n,  and  attrition  rate  p, 
the  number  of  attritions  y  is  described  by  the  Binomial 
(n,p)  distribution.  The  variance  of  the  estimator  of 
£  =  y/n  is  given  by  the  familiar  formula 

var  (£)  ■Ej-1-~JL)  .. 

n  (2.1) 
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The  Freeman-Tukey  double  arcsine  transformation  was  used 
by  Tucker  [Ref.  l:p.  55]  and  Robinson  [Ref.  2:pp.  74-75]  in 
order  to  map  the  raw  losses  to  the  transformed  scale.  This 
formula  is  given  by 


x  =  0.5[n  +  0 . 5] 1 >  2  I  sin- 1  [ 


2v_ 


(n  +  1) 


-  1] 


(2.2) 


+  sin“I  -  ui 


Note  that  n=n(t)  may  change  with  time  and  thus,  so  would  x. 
The  question  of  averaging  over  time  before  or  after  applying 
the  Freeman-Tukey  transformation  is  not  discussed  in  this 
paper.  One  of  the  goals  of  using  this  transformation  is  to 
stabilize  the  variance  of  this  estimate.  In  the  work  of 
Tucker  and  Robinson,  the  assumption  of  normally  distributed 
random  variables  with  common  variance  was  the  setting  in 
which  the  James-Stein  estimator  was  expected  to  perform 
well.  Indeed,  for  n  moderate  or  large  and  p  not  too 
extreme,  the  variance  of  the  Freeman-Tukey  transformation  is 
approximately  one.  But  as  discovered  by  Robinson,  this 
assumption  does  not  hold  well  when  applied  to  cells  with 
small  inventory.  A  more  careful  look  at  the  varience  of  the 
transformed  data  is  taken  here,  with  particular  interest  in 
the  region  of  unstable  variance  where  n  and  p  are  small. 

Some  exploratory  graphical  work  has  revealed  that  the 
natural  log  of  the  variance  of  x,  the  transformed  value. 
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behaves  quite  well  as  a  linear  function  of  u(  =E[X])  and  p- 
1  in  this  unstable  zone.  This  leads  to  the  interpolatory 
formula : 


var(x)  *  1.6835U-- •***  (p  -  l)-*8*i,  (2.3) 
1.001  S  p  S  2.2 

The  upper  limit,  p  =  2.2,  is  the  value  for  which  var(x)  =  1. 
Since  the  least  squares  fitting  process  was  applied  to  the 
logarithms,  p  can  not  be  allowed  to  fall  to  one  or  below. 
The  lower  limit,  p  -  1.001,  is  an  arbitrary  value  that  meets 
this  constraint.  The  coefficients  were  validated  by  a 
linear  regression  on  the  natural  log  of  this  formula.  The 
actual  variance  function  for  n  =  1,2,3,4,5,7,10  appear  in 
Figure  1,  and  the  fit  is  remarkably  good  for  n  i.  3.  The 
upper  graph  shows  the  region  of  special  interest  for  this 
paper  and  this  portion  of  the  curve  is  fitted  by  the  formula 
in  Equation  2.3  above.  The  lower  plot  exhibits  the  outer 
tale  effect  yet  to  be  explored.  The  small  values  of  n(<3) 
do  not  maintain  the  stable  variance  region,  i.e.,  var[x]=l, 
for  long  before  falling  back.  As  the  value  for  n  increases, 
the  more  stable  the  variance  becomes.  Table  2  is  a  partial 
display  of  the  residuals  for  selected  n  values,  fifteen 
equally  spaced  values  in  the  range  1.001  £  u  5  2.2,  computed 
by  the  linear  regression  formula  above.  For  values  of  n  >  2 
the  difference  between  actual  and  fitted  values  are 
considered  acceptable  for  this  work. 


9 


•  •  •  •  • 


Figure  1.  Variance  of  Freeman-Tukey  Transform  as 
Function  of  It's  Mean  for  Selected  n. 


Table  2.  RESIDUAL  LISTING  FOR  SELECTED  VALUES  OF  N 


N-I 

N-2 

N-4 

MM 

N  =  10 

0.07191 

0.03912 

0.02139 

0.017 

0.008313 

0.08466 

0.02569 

0.002853 

-0.001695 

-0.01106 

0.1033 

0.0239 

-0.005448 

-0.01049 

-0.0206 

0.122 

0.03002 

-0.00655 

-0.01185 

-0.02257 

0.1385 

0.04068 

-0.004026 

-0.009566 

-0.02067 

0.15 

0.05388 

0.0000676 

-0.005727 

-0.0166 

0.1552 

0.06853 

0.004027 

-0.002103 

-0.01291 

0.1527 

0.08234 

0.007447 

0.0009541 

-0.009981 

0.1414 

0.09481 

0.00948 

0.002557 

-0.008299 

0.1209 

0.1055 

0.009825 

0.002448 

-0.008128 

0.08958 

0.1129 

0.008369 

0.0004583 

-0.009929 

0.04832 

0.1169 

0.004997 

-0.003511 

-0.01327 

-0.005254 

0.1171 

0.0001116 

-O.OOS93 

-0.0182 

-0.06914 

0.113 

-0.006798 

-0.01654 

-0.02463 

-0.1466 

0.104 

-0.01534 

-0.02585 

-0.03366 

The  boundaries  for  u,  in  the  unstable  region,  as  a 
function  of  n  and  p  are  provided  in  Figure  2.  The  formula 
derived  was  quite  adequate  for  this  study. 

C.  ALTERNATE  TRANSFORM  INVERSION  METHODS 

There  are  three  issues  in  the  inversion  of  the  data  from 
the  transformed  space.  The  first  is  that  of  selecting 
values  of  n,  i.e.,  an  average  of  the  n(t)'s  appearing  in  the 
estimation  year  inventories.  This  issue  is  discussed  later 
in  the  chapter,  with  the  development  of  a  simulation  for 
comparing  mean  inventory  values.  The  second  is  that  of  the 
inversion  of  the  sum  of  the  two  arcsine  functions.  The 
third  involves  the  case  of  time  averaging  in  the  transformed 
scale  and  the  question  of  what  matching  form  of  average 
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INVENTORY 


12 


Figure  2.  Envelope  for  the  Non-Constant  Portion  of  the 

Fitted  Variance  as  Functions  of  Attrition  Rate 


inventory  should  be  used  in  the  inversion.  The  third  issue 
is  not  covered  in  this  paper,  and  is  still  pending. 
Treatment  of  the  second  issue  will  be  provided  first. 

The  method  used  in  the  pilot  studies  of  Tucker  and 
Robinson,  and  the  most  transparent  transform  inversion 
formula  is: 

*  0 . 5  [  1  -  sin  ( (n  +  *2.4) 

This  will  be  referred  to  as  the  Basic  Inversion.  The 
formula  offered  by  Carter  and  Rolph  [Ref.  4]  is  given  by: 


sin2  (0i )  + 


(1  -  2sin2  Qi ) 
( 4 (n/k)  +  2) 


!  Bi 


(2.5) 


where  0i  is  x/ (n  +  1/2) 1/2  and  Bi  is  the  estimated  amount  of 
shrinkage  in  the  transformed  scale  depending  on  the  Bayes 
model  used  (constant  or  proportional  prior) ,  k  is  the  number 
of  cells  and  n  is  inventory. 

Apparently,  Carter  and  Rolph  invert  prior  to  shrinking 
and  the  above  formula  provides  a  shrinkage  in  the  original 
scale.  Rewriting  the  above,  it  appears  to  shrink  towards 
£  ■  1/2,  and  is  given  by: 

-  1/2(4(n/io'v  2>  +  sin1  (Qi  )  (1  -4(n/**‘+  2>  (2.6) 

An  exact  Freeman-Tukey  inversion  formula  developed  by 
Miller  [Ref.  5]  was  recently  uncovered.  It  is  exact  in  the 
sense  that  for  a  single  cell,  the  formula  given  by: 
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p ( t'  =  0.5fl  -  sgn(cos  t) [1  -  (sin  t  + 

— I — I — |  (2.7) 

n 

where  t  »  x/(n  +  1/2) 1/2  +  n/2  and  a  is  the  average  central 
inventory,  returns  the  empirical  rate  y/n.  This  will  be 
referred  to  as  the  Freeman-Tukey  Exact  (FTE)  inversion. 
Since  this  formula  is  applied  using  an  average  value  of  n, 
there  may  be  concern  that  it  oscillates  between  successive 
integer  values  of  n.  Figure  3  indicates,  using  values  of 
n=5  and  10,  that  it  is  quite  smooth  as  an  interpolation 
formula.  In  the  article,  Miller  suggested  using  the 
harmonic  mean  of  the  n  values  provided  by  the  time  changes . 
This  issue  is  examined  later  in  this  chapter.  The  Miller 
formula  was  incorporated  in  the  Marine  Corps  Inventory 
Attrition  Analysis  Program  provided  in  Appendix  B. 

D.  SHRINKAGE  VISUALIZATION 

When  shrinkage  is  a  fixed  amount  for  all  cells,  it  is 
readily  visualized,  for  all  degrees  of  shrinkage,  by  the 
linear  diagrams  presented  in  the  works  of  Tucker  [Ref.  l:p. 
52]  and  Casella  [Ref.  6].  Robinson  [Ref.  2:p.  20]  also 
presented  diagrams  of  the  limited  translation  option  and  how 
the  shrinkage  is  affected.  All  of  these  diagrams  refer  to 
the  transformed  scale. 

Of  more  immediate  interest  is  the  question  of  how  to 
visualize  shrinkage  diagrams  of  this  type  in  the  original 
scale.  Illustrations  of  diagrams  of  this  type  are  provided 
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using  both  the  basic  inversion  formula  used  by  both  Tucker 


and  Robinson, 

and 

the 

FTE  formula. 

Figur 

e  4 

shows  how  the 

diagram  by  Casella 

[Ref 

6 :  p . 

2  4] 

looks 

in 

the  original 

scale  using 

the 

same 

value 

of 

n  for 

all 

cells.  An 

interesting  note  to  consider  is  the  effect  of  the  difference 
between  inversion  formulas,  i.e.,  in  the  graph  which 
combines  the  Basic  and  FTE  inverses,  the  FTE  inversion  tends 
to  shrink  at  a  slower  rate,  as  values  get  further  from  the 
mean.  For  a  modest  transition  into  the  case  of  non  constant 
n,  Figure  5  provides  for  two  values  of  n.  Finally,  Figure  5 
illustrates  the  point  with  many  different  values  of  n,  taken 
from  actual  Marine  Corps  officer  inventory  data.  In  this 
case,  note  that  Basic  inversion  shrinks  to  one  original 
scale  value  while  the  FTE  inverse  is  not  as  focused. 

E.  CHOICE  OF  AVERAGE  INVENTORY  VALUES  OVER  TIME 

A  Fortran  simulation  was  developed,  included  as  Appendix 
C,  to  test  the  performance  of  each  of  several  common  means 
(arithmetic,  geometric  and  harmonic)  in  the  comparison  of 
the  inverted  estimates  with  the  empirical  y/n.  Values  of  n 
were  produced  from  a  Poisson  distribution  with  rate  X  (1  S  X 
<  20),  and  attrition  values,  y,  were  produced  from  the 
Binomial  (n,p)  distribution  where  a  was  in  the  range  .01  s  p 
2:.  40.  The  n  and  y  values  were  then  transformed  by  using  the 
Freeman-Tuksy  formula.  An  exhaustive  study  of  inversion 
variations  with  both  the  FTE  formula  and  the  basic  formula 
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BASIC  NVERSC  (N-IO) 


Figure  4.  Non  Linear  Shrinkage  Diag 
Value  of  Inventory. 


BASIC  INVERSE  (N- 10.20) 


Figure  5.  Non  Linear  Shrinkage  Diagrams  for  Two  values 
of  Inventory. 


Figure  6.  Non  Linear  Shrinkage  Diagrams  with  Empirical 
Bayes  Estimators  for  Test  Case. 


used  by  Robinson  were  compared.  Table  3  contains  selected 
outputs  which  are  typical  of  the  broader  study  and  shows  how 
the  means  performed  against  each  other  in  three  different 
computations  listed  below. 

1.  The  BASIC  Method  utilizes  the  Basic  inversion  formula 
and  the  average  transformed  scale  value. 

2.  The  FTE1  Method  utilizes  Equation  2.7  with  the  t 

values  equal  to  the  average  (over  replications,  i.e., 
time)  of  x/(n  +  1/2) 1/2  and  x  given  by  Equation  2.2. 

More  explicitly,  see  Equation  A. 16. 

3.  The  FTE2  Method  utilize  Equation  2.7  with  the  t 
values  equal  to  the  average  (over  replication,  i.e., 
time)  of  x  directly,  and  x  given  by  Equation  2.2 
followed  with  a  division  by  (n  +  l/2)l/z 

In  the  context  of  this  study,  the  arithmetic  mean  tended 
to  perform  as  well  or  better  than  the  geometric  or  harmonic 
means.  The  choice  of  arithmetic  mean  also  circumvents  the 
problem  of  an  undefined  term  when  m  *  0  is  encountered 

while  using  the  harmonic  of  geometric  means. 


F.  QUESTIONS  RELATING  TO  VALIDATION 

The  pilot  studies  by  Tucker  and  Robinson  both  used  four 
years  (1977-80)  for  estimation  of  the  attrition  rates  and 
the  remaining  three  years  (1981-83)  for  validation.  The 
estimation  rates  were  used  in  an  attempt  to  forecast  the 
values  for  the  three  validation  years .  The  quality  of  the 
forecasts  deteriorates  as  lead  time  increases  in  the  future. 
This  was  anticipated,  in  fact,  Rowe  et  al . ,  at  NPRDC  have 
shown  a  time  series  effect  in  their  work  forecasting 
attritions  in  the  Navy  [Ref.  7]  and  are  currently 
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Table  3.  SIMULATION  OUTPUT  OF  MEAN  VALUES  FOR  INVERSION 
UTILIZATION  _ 


REPS 

P 

METHOD 

ARITHMETIC 

GEOMETRIC 

HARMONIC 

2.0 

40 

0.05 

BASIC 

0.1068 

0.1176 

0.1308 

FTE1 

0.0252 

0.01S9 

0.0122 

FTE2 

0.0181 

0.0208 

0.0245 

2.0 

40 

0.10 

BASIC 

0.1399 

0.1537 

0.1708 

FTE1 

0.0535 

0.0458. 

0.0369 

FTE2 

•  0.0489 

0.0558 

0.0648 

2.0 

40 

0.20 

BASIC 

0.2402 

0.2629 

0.2906 

FTE1 

0.1900 

0.1829 

0.1739 

FTE2 

0.1657 

0.1869 

0.2144 

2.0 

40 

0.40 

BASIC 

0.3713 

0.4038 

0.4429 

FTE1 

0.3442 

0.3403 

0.3353 

FTE2 

0.3325 

0.3715 

0.4213 

4.0 

40 

0.05 

BASIC 

0.0750 

0.0815 

0.0921 

FTE1 

0.0235 

0.0200 

0.0147 

FTE2 

0.0200 

0.0221 

0.0258 

4.0 

40 

0.10 

BASIC 

0.1280 

0.1387 

0.1564 

FTE1 

0.0768 

0.0724 

0.0655 

FTE2 

0.0726 

0.0798 

0.0922 

4.0 

40 

0.20 

BASIC 

0.1965 

0.2125 

0.2386 

FTE1 

0.1446 

0.1406 

0.1339 

FTE2 

0.1485 

0.1628 

0.1869 

4.0 

40 

0.40 

BASIC 

0.4382 

0.4694 

0.5185 

FTE1 

0.4564 

0.4558 

0.4549 

FTE2 

0.4280 

0.4638 

0.5223 

8.0 

40 

0.05 

BASIC 

0.0552 

0.05S0 

0.0616 

FTE1 

0.0181 

0.0166 

0.0148 

FTE2 

0.0156 

0.0165 

0.0177 

8.0 

40 

0.10 

BASIC 

0.1226 

0.1285 

0.1363 

FTE1 

0.0922 

0.0904 

0.0880 

FTE2 

0.0828 

0.0873 

0.0933 

8.0 

40 

0.20 

BASIC 

0.2035 

0.2131 

0.2255 

FTE1 

0.1778 

0.1762 

0.1741 

FTE2 

0.1708 

0.1797 

0.1915 

8.0 

40 

0.40 

BASIC 

0.3799 

0.3962 

0.4170 

FTE1 

0.3700 

0.3693 

0.3684 

FTE2 

0.3662 

0.3837 

0.4065 

recommending  a  second  order  auto  regression  scheme  for  this 
purpose.  In  this  regard,  the  work  of  Tucker  and  Robinson 
with  validation  may  be  viewed  as  forecasting  by  persistence , 
and  in  light  of  the  NPRDC  work,  can  hardly  be  expected  to 
perform  well  in  an  absolute  sense.  Their  work  does, 
however,  provide  some  comparisons  that  can  be  useful  in 
choosing  among  alternative  techniques ,  and  much  can  be 
learned  from  them. 

The  first  measure  considered  is  the  Mean  Squared  Error 
(MSE)  of  forecasts.  The  shrinkage  estimation  theory 
provides  that  these  values  should  be  about  unity  if 
persistence,  independence  over  time,  was  agreeable.  That 
is,  useful  validations  can  be  obtained  without  modeling  the 
time  dependence  of  the  attrition  process.  Only  for  a  very 
few  cases  does  this  hold  (i.e.,  Warrant  Officers,  see 
Tucker) .  However,  if  one  method  of  estimation  consistently 
produces  noticeably  smaller  MSE  values  than  do  the  other 
methods,  this  would  provide  sufficient  support  for  that 
method's  continued  development. 

The  second  measure  is  designed  to  consider  performance 
in  the  original  scale  and  is  patterned  after  a  chi-square 
statistic 


[actuals  -  forecasts)2 
forecasts 


If  the  k  cells  are  independent,  then  equation  2.8  is 
approximately  a  chi-square  random  variable  with  k  degrees  of 
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freedom.  Hence,  the  man  value  Is  about  k  and  the  standard 
deviation  is  (2k)1/2.  Once  again,  we  generally  do  not  see 

i 

values  this  small,  but  use  it  for  comparison  purposes.  As  a 
side  note,  the  value  of  k  when  using  the  MLE  is  generally 
smaller  than  that  of  the  other  estimators  due  to  the 
omission  of  those  cells  with  MLE  equal  to  zero  over  all 
estimation  years,  i.e.,  zero  forecasts.  This  provides  and 
unnatural  advantage  for  the  MLE,  which  must  be  kept  in  mind 
when  Interpreting  the  results. 

Two  disturbing  features  have  emerged  from  the 
exploratory  analysis.  Often  a  technique  that  performs  well 
for  some  cases  using  the  MSE  measure  is  not  comparably 
supported  by  the  chi-square  measure.  It  seems  the  two 
measures  do  not  weight  the  common  features  equally,  and 
further  study  is  necessary  to  develop  understanding  of  the 
results . 

Discussion  with  personnel  at  NPRDC  concerning  these 
matters  has  led  to  the  introduction  of  a  third  measure, 
which  will  be  called  the  Mean  Absolute  Deviation  (MAD). 
This  value  is  computed  in  the  original  scale.  Moreover, 
since  the  cost  structure  of  over-forecasting  is  not 
symmetric,  we  separate  the  two  parts  into  the  average 
overage  deviation  given  by: 

•t 

1  k 

ODEV  ■  £  71  [forecasts  -  actuals] ♦  (2.9) 

k  X 
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and  average  underage  deviation 


UDEV  -  *■  E  [actuals  -  forecasts]*  (2.10) 

K  1 

where  the  "+"  refers  to  the  positive  part,  as  computed 
separately.  Finally, 

MAD  -  ODEV  +  UDEV  (2.11) 

This  measure  also  has  shown  that  in  most  cases,  the 
overage  exceeds  the  underage  value,  which  in  turn  may  lead 
to  forecasting  the  need  for  an  excess  number  of  Marine 


officers . 


III. 


EMPIRICAL  BAYES  ESTIMATION 


A.  BACKGROUND 

Although  the  studies  of  Tucker  and  Robinson  showed 
improvement  over  the  current  methods,  no  clear  procedure  was 
established  and  no  plans  for  implementation  were  drafted. 
Problem  areas  were  identified  while  working  with  small  cell 
inventories  and  low  loss  rates.  In  the  transformed  scale, 
the  James-Stein  assumption  of  normal  distribution  of  the 
observations  with  equal  variance  did  not  hold  well  in  such 
cases.  An  alternate  method  of  loss  rate  estimation  is 
required  to  deal  with  this  special  case. 

B .  INTRODUCTION 

Empirical  Bayes  is  simply  a  generalization  of  the  James- 
Stein  estimator  for  use  with  small  cell  inventories  and 
unstable  variance.  Random  variables  are  assumed  to  be  of 
the  form: 


Xi  '  N(9i ,di ) ,  i-1,2 _ _ k  (3.1) 

where  0i  are  the  unknown  parameters  to  be  estimated  and  di 
are  known  variances.  The  Bayesian  assumption  is  that  the 
©t  's  are  themselves  from  some  distribution,  in  this  case 
also  from  the  normal  distribution  with  mean,  v,  and 

proportionality  constant  pi  is 


variance,  pi  a..  The 
normalized  so  that 
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S  p  i  =  k , 


■i  —  1  O 

—  "*  f  «  f  •  • 


(3.2) 


•  ,k 

where  k  is  the  number  of  observations  or  cells.  Thus, 

9i  "  Nfv.pia),  i=l,2 . k.  (3.3) 

This  distribution  is  known  as  the  prior  distribution  of  ©t 
and  represents  the  additional  assumed  information  on  the 
value  of  ©  before  observing  X  whose  distribution  depends  on 
©  [Ref.  8:p.  553].  The  use  of  this  information  produces  an 
estimate  wh: is  a  weighted  average  of  the  prior  mean,  v, 
and  the  sample  estimate  x.  In  the  empirical  Bayes  context, 
the  values  of  v  and  a,  are  also  estimated  from  the  data 
since  all  the  information  necessary  is  contained  in  the 
marginal  distribution  of  Xi  (unconditional  on  9i )  [Ref.  6:p. 
83] .  This  chapter  focuses  on  that  method. 

For  a  more  detailed  explanation  of  the  Bayesian  family 
of  estimators,  the  reader  is  referred  to  Casella  [Ref.  6], 
Carter  and  Rolph  [Ref.  4]  and  Efron  and  Morris  [Ref.  9]. 

C.  THE  ESTIMATION  METHOD 

Fay  and  Herriot  [Ref.  10]  discuss  the  empirical  Bayes 
technique  and  provide  an  estimate  for  the  prior  mean  as 

S  xi  /  (a*  +  di  ) 

V*  «  _J _  ,  i=l ,  2 . k  (3.4) 

z:  1/ (a*  +  di  ) 

where  a*  is  an  unbiased  estimate  of  a.  The  empirical  Bayes 
estimate  of  0,  commonly  called  the  mean  of  the  posterior 
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where  a 


is  an  unbiased  estimate  of  a.  The  empirical  Bayes 
estimate  of  0,  commonly  called  the  mean  of  the  posterior 
distribution  becomes 


5E® 


( a* 


+  di  ) 


)  Xl 


dt 


(a*  +  di  ) 


)  v* 


i-1,2 . k  (3.5) 


This  is  the  mean  of  the  conditional  distribution  of  0  given 
x  [Ref.  8:p.  556].  The  weights  of  each  of  the  prior  and 
sample  means  are  determined  by  the  prior  and  observed 
variances . 

As  a  measure  of  the  worth  of  this  estimator,  when  che 
number  of  x's  is  £  4  [Ref.  11],  the  risk,  or  sum  of  squared 
differences  between  the  true  unknown  parameter  and  the 
empirical  Bayes  estimate,  is  less  than  that  using  the 
observed  x’s,  for  all  8i .  Efron  and  Morris  [Ref.  9]  provide 
a  rigorous  proof. 

An  important  result  of  the  empirical  Bayes  estimator  is 
that  shrinkage  values  are  no  longer  constant,  but  depend 
upon  information  contained  in  each  cell.  Additionally,  the 
point  toward  which  the  cell  values  are  shrunk  is  the 
weighted  mean,  not  the  grand  mean  of  the  observations. 
Appendix  A  provides  the  derivation  of  the  empirical  Bayes 
estimator,  variable  shrinkage  factors  and  calculation  of 
the  prior  variance  used  in  this  study  as  they  pertain  to  the 
Marine  Corps  officer  attrition  rates. 
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IV.  VALIDATION  PROCEDURES  AND  RESULTS 


A .  GENERAL 

The  outputs  displayed  in  this  chapter  are  excerpts  from 
that  produced  by  the  Marine  Corps  Personnel  Inventory 
Attrition  Analysis  (MCIAAP)  (see  Appendix  B) .  Initially, 
comparisons  will  be  made  with  the  output  displayed  by 
Robinson  to  identify  the  validity  of  the  Fortran  code  with 
the  APL  code  used  by  Robinson,  then  new  approaches  will  be 
displayed  in  an  attempt  to  improve  the  results.  The  Fortran 
code  in  the  MCIAAP  incorporates  most  of  the  procedures 
studied  by  Tucker  and  Robinson.  Specifically,  the  MLE , 
Transformed  Scale  Cell  Average  (TSCA)  and  the  James-Stein 
(JS)  estimators,  shrinkage  procedures  and  a  basic  inversion 
formula.  The  FTE  inversion  formula,  the  new  measures  of 
effectiveness  from  NPRDC ,  are  also  included.  For  the  MAD 
measure  of  effectiveness,  only  maximum  likelihood  and  the 
modified  empirical  Bayes  variables  will  be  presented,  since 
these  showed  the  most  promise  over  all.  All  other  variables 
are  displayed  by  the  MCIAAP  if  interest  is  warranted  by  the 
reader.  Chapter  II  describes  all  of  these  methods  in  some 
detail . 

B.  HISTORICAL  COMPARISONS 

Taoles  4-9  were  developed  using  the  same  aggregation 
scheme  and  format  as  Robinson  [Ref.  2:pp.  39-41],  although 
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not  all  of  the  variables  are  included  (MLE,  TSCA,  and  James 
Stein  are  the  only  variables  duplicated  in  this  paper) . 
Additionally,  Robinson  presented  the  number  of  cells,  as  k 
values,  used  for  evaluation  of  each  estimate  in  the 
aggregate  Figure  of  Merit  (FOM)  display.  These  k  values  are 
not  repeated  in  the  current  FOM  tables,  but  are  included,  as 
parenthetical  values,  and  discussed  in  the  new  MAD  measure 
of  effectiveness  tables.  Since  this  is  the  case  here,  the 
reader  is  reminded  that  the  apparent  competitive  performance 
of  the  MLE  in  the  following  cases  is  due,  in  part,  to  the 
possible  advantage  of  having  a  smaller  number  of  cells  used 
in  this  MOE's  computation. 

Since  the  computations  were  identical  for  both  studies 
for  these  tables,  the  output  was  expected  to  be  nearly 
identical  or  some  pattern  discernable  between  the  two.  In 
very  few  cases  was  this  found  to  be  true,  and  after 
verifying  that  the  algorithms  coded  both  in  the  APL  and 
Fortran  versions  were  correct,  we  have  no  reasonable 
explanation  for  the  discrepancies  of  the  programs.  It  is 
suspected  that  the  transfer  of  raw  data  into  an  APL 
workspace,  accomplished  by  Tucker  and  Robinson,  may  be  at 
variance  with  that  into  Fortran  files,  done  by  contractor. 

1.  Aviation  Aggregate  Case 

In  the  case  of  aviation  FOM ' s ,  displayed  in  Table  4, 
1981  values  for  LtCol's  in  the  transformed  scale  are 
identical  to  Robinson’s  [Ref.  2:p.  39]  while  the  last  two 
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Table  4.  AVIATION  FIGURES  OF  MERIT 


19S1  1982  1983 

TRANSFORMED  FOM 


1st  Lt 


MLE 

4.0863 

10.6590 

11.1338 

TSCA 

3.7420 

10.3636 

10.8711 

JS 

3.9694 

10.5671 

11.1626 

MOD  JS 

4.7352 

11.2615 

,  12.1206 

EB 

4.0753 

10.7015 

11.3303 

MOD  EB 

3.1833 

5.0568 

6.2870 

Lt  Col 

MLE 

4.3664 

8.5033 

8.1506 

TSCA 

5.7768 

10.2994 

10.3353 

JS 

5.  /3b4 

10.2468 

10.2972 

MOD  JS 

5.7208 

10.1865 

10.2756 

EB 

5.8482 

10.2298 

10.3534 

MOD  EB 

4.3301 

3.9824 

2.7186 

ORIGINAL  FOM 


1st  Lt 


MLE 

21.0501 

43.1546 

53.0538 

TSCA 

56.4603 

93.7810 

98.4715 

TSCAM 

53.7S84 

90.5164 

95.3329 

JS 

58.0150 

95.2150 

100.1664 

JSM 

55.3844 

91.9406 

97.0163 

MOD  JS 

63.0345 

99.S419 

105.4S97 

MOD  JSM 

60.6413 

96.5615 

102.3172 

EB 

58.8824 

96.2449 

101.24S2 

EBM 

56.2111 

92.9699 

98.0958 

MOD  EB 

19.7683 

38.4637 

49.7393 

MOD  EBM 

20.7987 

36.3673 

48.1659 

Lt  Col 

MLE 

33.9315 

59.4200 

22.1156 

TSCA 

38.4478 

57.0987 

69.0641 

TSCAM 

35.0390 

48.9374 

59.6957 

JS 

37.8248 

55.1913 

67.1696 

JSM 

33.2016 

46.9901 

56.7135 

MOD  JS 

38.0291 

51.5325 

63.3908 

MOD  JSM 

34.5466 

45.1554 

52.8285 

EB 

40.3421 

54.9389 

67.7091 

EBM 

36.7498 

46.9657 

57.8567 

MOD  EB 

40.2056 

70.5381 

60.5749 

MOD  EBM 

185.7476 

65.1553 

168.3976 

years  are  smaller.  The  IstLt  FOM's  are  larger  for  all  three 
years.  The  modified  empirical  Bayes  (Mod  EB!  values  were 
best  overall  yet  still  may  be  too  large  (ideal  values  tend 
to  be  close  to  unity) .  It  is  interesting  to  note  that  in 
most  cases,  the  risk  values  tend  to  increase  with  future 
forecasts,  yet  the  Mod  EB  FOM  decreases  as  future  forecasts 
increase  with  time. 

In  the  original  scale,  MLE  tended  to  be  smaller  chan 
Robinson's  for  all  years,  while  TSCA  and  James-Stein  (JS) 
were  all  larger.  The  FTE  inversion  worked  best  for  all 
variables  except  for  Mod  EB  in  the  first  and  third 
validation  years.  MLE  and  Mod  EB  worked  best  for  IstLt 's 
and  are  approximately  within  the  range  of  one  standard 
deviation  of  the  number  of  cells  in  the  computation.  ^or 
LtCol's,  FTE  inverted  JS  was  best  with  MLE  a  competitive 
second,  yet  all  were  not  in  the  desired  range. 

In  the  MAD  measure,  Table  5,  results  are  mixed 
between  future  forecasts  (no  set  pattern)  and  variables  for 
the  two  ranks  (FTE  Mod  EB  was  best  for  IstLt 's  and  MLE  was 
best  LtCol's).  The  mix  in  overage  and  underage  is  most 
desirable  in  the  first  year  for  all  cases,  yet  the  number  of 
cells  (k)  forecasted  as  overage  exceed  the  underage  forecast 
cell  numbers  in  most  cases.  This  is  cause  for  concern  and 
the  performance  of  these  estimates  in  this  aggregate  is 
suspect.  The  second  year  values  for  LtCol  estimates 
displays  a  preferred  ratio  of  k  values,  but  the  total  MOE 
value  may  still  be  too  large. 
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Table  5.  MEAN  ABSOLUTE  DEVIATION  (AVIATION) 


1981 

1982 

1983 

1st  Lt 

MLE 

UNDERAGE 

0.1256C  3) 

0.0439(  1) 

0.3069(  1) 

OVERAGE 

2.3985(  8) 

3.9267(10) 

4.2522(10) 

TOTAL 

2.5240(11) 

3.9706(11) 

4.5591(11) 

MOD  EB 

< 

UNDERAGE 

0. 1 822(  4) 

0.0410(  I) 

0.32S  I(  1) 

OVERAGE 

2.3625(  8) 

3.9040(11) 

4.2351(11) 

TOTAL 

2.5447(12) 

3.9451(12) 

4.5633(12) 

MOD  EBM 

UNDERAGE 

0.25 1S(  5) 

0.0580{  2) 

0.3420(  I) 

OVERAGE 

2.1654(  7) 

3.6427(10) 

3.9677(11) 

TOTAL 

2.4172(12) 

3.7007(12) 

4.3098(12) 

Lt  Col 

MLE 

UNDERAGE 

0.6005(  5) 

1.5139(10) 

0.0396(  1) 

OVERAGE 

0.7925(  8) 

0.3928(  3) 

1.3467(12) 

TOTAL 

1.3930(13) 

1.9067(13) 

1.3863(13) 

MOD  EB 

UNDERAGE 

0.6053(  6) 

1.4717(  9) 

0.084 1(  2) 

OVER-AGE 

1.3083(10) 

1.0970(  8) 

2.2331(14) 

TOTAL 

1.9136(16) 

2.56S7(17) 

2.3173(16) 

MOD  EBM 

UNDERAGE 

0.6S90(  6) 

1.5366(10) 

0.1 167(  2) 

OVERAGE 

1.0120(10) 

0.71 79(  7) 

1.7501(14) 

TOTAL 

1.7010(16) 

2.2546(17) 

l.S669(  16) 

2.  Combat  Support  Aggregate  Case 

Table  6  displays  the  combat  support  aggregate  output 
from  the  program.  In  the  transformed  scale  comparison  with 
Robinson  [Ref.  2:p.  401,  the  values  for  the  two  ranks  take 
opposite  directions.  Where  the  IstLt  FOM  value*  are  all 
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Table  6.  COMBAT  SUPPORT  FIGURES  OF  MERIT 


TRANSFORMED  FOM 


1st  Lt 


MLE 

2.1877 

2.9028 

2.2447 

TSCA 

1.6329 

2.4580 

1.8696 

JS 

1.6925 

2.4816 

1.8987 

MOD  JS 

2.3684 

3.0520 

2.5356 

EB 

2.3323 

3.1991 

2.6S68 

MOD  EB 

2.2293 

2.0048 

'  2.1531 

iol 

MLE 

1.3103 

1.6978 

2.3108 

TSCA 

0.7502 

1.4070 

1.7540 

JS 

0.6421 

1.2345 

1.5429 

MOD  JS 

0.7009 

1.0929 

1.3082 

EB 

0.7736 

1.2589 

1.5799 

MOD  EB 

1.8691 

2.1236 

4.9017 

MLE 

ORIGINAL  FOM 

111.4260 

112.1917 

70.2767 

TSCA 

73.5008 

90.0523 

73.9380 

TSCAM 

79.0947 

93.8244 

80.9764 

JS 

76.9945 

89.S064 

71.5890 

JSM 

429.8503 

84.4175 

65.2639 

MOD  JS 

143.1700 

113.0083 

90.6783 

MOD  JSM 

168.9385 

114.7083 

92.6599 

EB 

102.9525 

119.0900 

96.3373 

EBM 

108.3578 

111.2701 

86.8138 

MOD  EB 

6S.7251 

63.9054 

4S.2719 

MOD  EBM 

176.3961 

59.5760 

43.0037 

!ol 

MLE 

36.0152 

66.0261 

42.4097 

TSCA 

61.7769 

52.0346 

99.7917 

TSCAM 

65.9535 

49.6307 

107.1016 

JS 

61.0511 

49.2280 

96.7861 

JSM 

66.3509 

47.8490 

118.2564 

MOD  JS 

66.5470 

49.9366 

99.2566 

MOD  JSM 

81.4722 

50.2637 

103.8460 

EB 

69.7508 

55.9812 

104.5311 

EBM 

71.5728 

50.0281 

109.2012 

MOD  EB 

62.6316 

50.3125 

96.8693 

MOD  EBM 

83.1337 

48.6576 

127.1541 

•< 
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larger  (except  for  MLE  in  1981), the  LtCol  values  are 
smaller.  Overall,  all  estimates  are  competitive  with  each 
other  over  all  validation  years  with  TSCA  performing  best 
for  the  IstLt's  and  James-Stein  outperforming  the  rest  in 
the  LtCol' s.  Note  here  also  that  the  Mod  EB  estimate  did 
not  fair  well  in  either  case,  and  was  the  worst  estimate  in 
risk  for  the  LtCol 's. 

In  the  original  scale  ,  a  different  and  more  random 
effect  is  noticed  between  Table  6  and  the  results  of 
Robinson.  Risk  values  tended  to  fluctuate  with  forecast 
lead  time  between  estimates  for  IstLt's  where  the  LtCol 
values  were  more  stable,  yet  were  quite  different  from  the 
transformed  scale  pattern.  In  this  case,  Mod  EB  performed 
well  for  IstLt's,  and  was  the  only  one  within  the  desired 
range  of  standard  deviations.  MLE  for  the  LtCol ' s  was  also 
the  only  estimate  desirable  for  this  measure.  The  Miller 
inverse  showed  poorly  in  this  aggregate  for  all  estimation 
methods,  in  every  case  exceeding  it's  counterpart. 

In  the  MAD  measure  of  effectiveness,  Table  7,  MLE 
performed  the  best  in  both  ranks.  The  Mod  EB  effect  is 
competitive  in  the  IstLt's  case,  and  all  variables  in  the 
LtCol  case  are  closely  grouped.  The  number  of  cells  are 
heavily  weighted  towards  an  overage  forecast  in  all  three 
validation  years,  throughout  this  aggregate.  Once  again, 
this  is  not  a  desirable  quality  for  a  MOE. 
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Table  7.  MEAN  ABSOLUTE  DEVIATION  (COMBAT  SUPPORT) 


1981 

1982 

1983 

1st  Lt 

MLE 

UNDERAGE 

0.3737(14) 

0.1617(8)  . 

0.178K  8)  • 

OVERAGE 

0.6731(19) 

0.9574(25) 

0.8835(25) 

TOTAL 

1.0468(33) 

1.1191(33) 

1.0616(33) 

MOD  EB 

UNDERAGE 

0.4028(15) 

0.1275(  9) 

0.1427(  9) 

OVERAGE 

0.7311(40) 

0.99^2(45) 

0.9422(38) 

TOTAL 

1.1339(55) 

1.1246(54) 

1.0849(47) 

MOD  EBM 

UNDERAGE 

0.4592(16) 

0.1603(  9) 

0.1856(10) 

OVERAGE 

0.5930(39) 

0.S3 18(45) 

0.7575(37) 

TOTAL 

1.0522(55) 

0.9921(54) 

0.9431(47) 

Lt  Col 

MLE 

UNDERAGE 

0. 1 360(  9) 

0.16S8(  7) 

0.054 1(  5) 

OVERAGE 

0.1510(14) 

0.2591(16) 

0.3415(18) 

TOTAL 

0.2S70(23) 

0.4280(23) 

0.3957(23) 

MOD  EB 

UNDERAGE 

0.1638(11) 

0.1721(  7) 

0.1 162(  6) 

OVERAGE 

0.2232(29) 

0.3606(35) 

0.454S(29) 

TOTAL 

0.3S70(40) 

0.5327(42) 

0.5710(35) 

MOD  EBM 

UNDER-AGE 

0.2044(11) 

0.1901(  7) 

0. 1 40S(  6) 

OVERAGE 

0.1263(29) 

0.2240(35) 

0.2589(29) 

TOTAL 

0.3307(40) 

0.4142(42) 

0.3996(35) 

3.  Ground  Combat  Aggregate  Case 

In  the  final  aggregate  studied  by  Robinson  [Ref . 
2:p.  41],  that  of  ground  combat  MOS ' s ,  the  results  again,  as 
displayed  in  Table  8,  were  mixed  between  computational 
methods.  In  the  original  scale,  value  differences  tended  to 
be  larger  than  Robinson's  in  all  cases  for  IstLt's,  showing 
a  proportional  growth  over  time.  Conversely,  LtCol  values 
were  smaller  in  the  first  two  validation  years  and  slightly 
larger  in  the  final  year.  TSCA  looks  good  for  the  IstLt's 
with  the  others  being  competitive,  yet  the  size  of  the 
values  may  not  be  desirable.  Mod  EB  turned  out  the  best  for 
LtCol ' s ,  but  all  other  estimates  are  competitive  and  in  the 
desired  range. 

On  the  original  scale,  all  values  for  the  estimates 
for  IstLt's  are  too  large,  with  Mod  EB  performing  well  in 
that  regard.  LtCol  values  are  all  larger  than  Robinson's, 
but  the  estimates  are  all  within  range  of  one  standard 
deviation  and,  therefore,  competitive.  MLE ,  once  again, 
shows  well  for  this  rank.  In  opposing  fashion  to  that  of 
the  combat  support  aggregate,  the  FTE  inversion  outperformed 
all  Basic  inversions  for  all  variables. 

The  MAD  measure  of  effectiveness.  Table  9  confirms 
the  results  found  above  in  that  Mod  EB  was  best  for  IstLt's 
and  MLE  did  well  for  LtCol' s.  The  values  for  IstLt  were 
found  to  be  too  large  overall,  while  all  estimates  for  LtCol 
were  competitive  and  in  the  desired  range.  The  mix  between 
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Table  8.  GROUND  COMBAT  FIGURES  OF  MERIT 


1981  1982  1983 

TRANSFORMED  FOM 


1st  Lt 


MLE 

2.5764 

4.4136 

6.S744 

TSCA 

1.9741 

3.8072 

6.1152 

JS 

2.1737 

3.9738 

6.3970 

MOD  JS 

3.0859 

4.7919 

7.5619 

EB 

2.5558 

4.3626 

.  6.9607 

MOD  EB 

3.1538 

3.5202 

4.7069 

Lt  Col 

MLE 

1.2126 

1.9674 

3.1084 

TSCA 

1.5715 

2.3553 

3.5566 

JS 

1.5571 

2.2719 

3.4365 

MOD  JS 

1.6227 

2.1317 

3.1S83 

EB 

1.7087 

2.2464 

3.3345 

MOD  EB 

1.1964 

2.5830 

2.5029 

ORIGINAL  FOM 


1st  Lt 


MLE 

87.3362 

119.2637 

143.9128 

TSCA 

87.9712 

116.6767 

145.6236 

TSCAM 

83.7648 

111.3995 

142.0902 

JS 

95.4717 

122.9425 

152.7354 

JSM 

92.0238 

117.5264 

149.3965 

MOD  JS 

124.6490 

147.6611 

180.5502 

MOD  JSM 

134.9472 

142.5461 

178.6173 

EB 

109.6349 

136.6321 

167.3234 

EBM 

105.5708 

131.0400 

163.8189 

MOD  EB 

64.3647 

81.2071 

101.9511 

MOD  EBM 

63.6623 

77.8159 

100.7150 

Lt  Col 

MLE 

35.0538 

45.8172 

53.0111 

TSCA 

55.9649 

198.8659 

123.1306 

TSCAM 

43.5234 

201.9953 

122.8648 

JS 

56.3346 

198.9037 

122.2786 

JSM 

43.9957 

201.3299 

115.7143 

MOD  JS 

60.0730 

202.2959 

122.9959 

MOD  JSM 

52.3940 

210.7551 

121.4254 

EB 

62.7441 

205.4249 

126.6544 

EBM 

50.4049 

209.3853 

125.7642 

MOD  EB 

48.0871 

197.7316 

111.9375 

MOD  EBM 

38.1231 

323.3654 

242.4398 

the  number  of 


over-estimates  and  under-estimates,  i.e., 
cells  forecasted  over  or  under  the  actuals,  is  not  a 


desirable 

one,  yet  the 

same 

mix 

holds  true 

for 

all 

validation 

years.  This 

also 

was 

the  finding 

of 

both 

preceding  aggregates. 


Table  9.  MEAN  ABSOLUTE  DEVIATION  (GROUND  COMBAT) 


1981 

1982 

1983 

1st  Lt 

MLE 

UNDERAGE 

0.1525(  5) 

0.1855(10) 

0.2466(10) 

OVERAGE 

1.9322(30) 

2.5979(25) 

3.1691(25) 

TOTAL 

2.0847(35) 

2.7834(35) 

3.4157(35) 

MOD  EB 

UNDERAGE 

0.2340(10) 

0.1323(  8) 

0.2220(  8) 

OVERAGE 

1.7743(33) 

2.3086(32) 

2.8449(30) 

TOTAL 

2.0083(43) 

2.4409(40) 

3.0669(38) 

MOD  EBM 

UNDERAGE 

0.2618(10) 

0.1705(10) 

0.2617(10) 

OVERAGE 

1.6492(33) 

2. 1 73S(  30) 

2.6935(2S) 

TOTAL 

1.9111(43) 

2.3443(40) 

2.9552(38) 

Lt  Col 

MLE 

UNDERAGE 

0.2313(10) 

0.3034(11) 

0. 1 36S( 10) 

OVERAGE 

0.4182(25) 

0.4202(24) 

0.6453(25) 

TOTAL 

0.6495(35) 

0.7236(35) 

0.7821(35) 

MOD  EB 

UNDERAGE 

0.2315(11) 

0.3299(14) 

0.1426(12) 

OVERAGE 

0.6392(34) 

0.6922(32) 

0.9523(31) 

TOTAL 

0.8706(45) 

1.0221(46) 

1.0950(43) 

MOD  EBM 

UNDERAGE 

0.2650(12) 

0.3790(15) 

0.1864(12) 

OVERAGE 

0.4071(33) 

0.4427(31) 

0.6535(31) 

TOTAL 

0.6720(45) 

0.S2 17(46) 

0.8399(43) 
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4.  Test  Aggregate  Case 

In  an  effort  to  experiment  with  aggregate  groupings, 
another  case  was  built  upon  which  all  estimation  methods 
were  applied,  and  the  results  compared  with  previously- 
studied  groups.  The  grade  was  held  fixed,  in  this  case  the 
grade  of  Captain  was  used,  with  a  LOS  in  the  eighth  year. 
Several  MOS ' s  were  chosen  so  that  to  insure  an  aggregate 
with  non-zero  inventories.  In  the  transformed  scale,  Table 
10,  all  estimates  have  small  and  competitive  values,  with 
Mod  EB  and  MLE  doing  well.  In  the  original  scale,  once 
again,  all  estimates  are  well  within  range,  but  Mod  EB  and 
MLE  were  the  standouts.  The  FTE  inversion  option  did  not 
fair  well  with  this  aggregate.  All  estimates  in  the  MAD 
measure  of  effectiveness,  Table  11,  are  very  competitive, 
and  the  first  validation  year  mix  (overage  and  underage)  may 
be  the  best  yet.  It  is  interesting  to  note  that  the  cell 
split  in  the  first  validation  year  between  overage  and 
underage  forecasts  is  contrary  to  that  of  the  historical 
displays  above.  Although  the  total  number  of  cells  used 
(k=6)  is  much  smaller  than  the  previous  table,  it  is  felt 
that  the  aggregate  scheme  was  the  major  factor  in  this 
finding,  and  the  results  are  promising. 
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Table  10.  TEST  AGGREGATE  FIGURES  OF  MERIT 


1981 

1982 

1983 

TRANSFORMED  FOM 

Captain 

• 

MLE 

0.9603 

2.3463 

2.0842 

TSCA 

1.4151 

3.3033 

2.9152 

JS 

1.3753 

3.3969 

2.8725 

MOD  JS 

1.3943 

3.8165 

2.8829 

EB 

1.6840 

4.2698 

3.2957 

MOD  EB 

0.8249 

1.7132 

1.2980 

ORIGINAL  FOM 

Captain 

MLE 

3.0402 

4.8104 

4.9399 

TSCA 

4.6161 

10.2779 

7.5714 

TSCAM 

8.1467 

7.6241 

8.8347 

JS 

4.3970 

10.3000 

7.4726 

JSM 

7.3199 

7.6666 

8.3860 

MOD  JS 

4.2990 

10.7640 

7.6371 

MOD  JSM 

6.0079 

8.3070 

7.9259 

EB 

5.3660 

12.0337 

8.8679 

EBM 

6.7181 

9.3915 

8.7959 

MOD  EB 

1.4463 

6.0323 

3.9066 

MOD  EBM 

4.6465 

4.3354 

5.9175 

Table  1 1.  MEAN  ABSOLUTE  DEVIATION  (TEST  CASE) 


1981 

1982 

Captain 

MLE 

UNDERAGE 

0.2625(  3) 

0.0463(  1) 

OVERAGE 

0. 1 35 1(  3) 

0.7374(  5) 

TOTAL 

0.3976(  6) 

0.7837(  6) 

MOD  EB 

UNDERAGE 

0.1342(  4) 

0.0000(  0) 

OVERAGE 

0.2575(  2) 

0.9601{  6) 

TOTAL 

0.39 18(  6) 

0.960 1(  6) 

MOD  EBM 

UNDERAGE 

0.3712C  4) 

0.0809(  1) 

OVERAGE 

0.0852(  2) 

0.6148(  5) 

TOTAL 

0.4563(  6) 

0.6957(  6) 

.  * 


1983 


0.2097(  2) 
0.439 1(  4) 
0.6488(  6) 


0.0965(  1) 
0.6306(  5) 
0.7270(  6) 


0.2352(  2) 
0.3025(  4) 
0.5377(  6) 
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v.  .  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  RESULTS 

Omitting  the  fact  of  differences  in  computational 
methods  between  this  paper  and  the  Robinson  thesis,  it  is 
clear  that  the  results  are  still  inconclusive  at  this  time. 
In  the  historical  comparisons,  the  estimates  studied  varied 
in  their  performance  from  one  aggregate  to  another.  The 
different  risk  calculations,  in  the  transformed  and  original 
scales,  as  well  as  the  different  measures  of  effectiveness 
showed  little  agreement  in  estimation  methods.  Long  range 
estimation  still  remains  a  problem  since  with  time,  the 
estimate  risks  increase  significantly.  Much  emphasis  was 
put  on  the  modified  empirical  Bayes  estimator  in  this  paper 
and  those  specific  results  look  promising  under  possibly 
different  and  more  stable  aggregate  conditions. 

One  area  of  mention  here  concerning  the  results  has  to 
do  with  the  composition  of  inventories  used  in  the  selected 
aggregates.  It  was  noticed  that  the  range  of  inventory 
(smallest  to  largest  value  in  any  particular  MOS)  was  quite 
large,  and  may  be  affecting  the  results.  In  fact,  Carter 
and  Rolph  (Ref.  4:p.  883]  managed  this  point  by  further 
disaggregating  into  levels  of  activity  as  it  concerned  their 
work  with  fire  alarm  estimation. 
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B.  CONCLUSIONS 


The  results  of  the  previous  feasibility  studies  and  this 
paper  leads  to  several  conclusions.  The  chosen  aggregate 
still  seems  to  be  the  key  element  in  the  work  of  attrition 
estimation.  It  does  seem  clear  that  the  aggregation  scheme 
should  be  focused  on  fixed  ranks,  in  narrow  LOS  ranges  and 
multiple  MOS ' s .  This  conclusion  is  based  upon  the  relative 
success  of  the  test  aggregate  utilized  in  this  work,  and 
upon  similar  success  by  Carter  and  Rolph.  Additionally,  the 
partitioning  of  inventory  levels  in  the  aggregate  should 
stabilize  the  estimation  method  performance.  Ail  the 
estimation  method  studied  to  date  have  their  strengths,  but 
still,  no  one  method  of  estimation  excels  over  the  others. 

C .  RECOMMENDATIONS 

No  estimation  method  studied  in  this  paper  is 
recommended  for  implementation  at  this  time.  In  addition, 
no  method  is  recommended  for  removal  from  consideration 
since  basic  aggregation  problems  are  felt  to  be  the  largest 
contributory  factor  to  the  mixed  results  found  here.  These 
recommendations  are  presented  for  further  study. 

1.  The  aggregation  method  offered  by  Captain  Larsen 
[Ref.  3]  should  be  examined  further  with  its 
relevance  to  the  scheme  in  the  test  case  used  here. 

2.  Utilization  of  the  detailed  data  provided  by  NPRDC  as 
soon  as  possible  in  order  to  expand  the  available 
information  and  open  new  avenues  of  aggregation. 

3.  Re-examine  the  inventory  mean  value  used  in  the  FTE 
inverse  to  determine  its  future  role  in  the 
estimation  process. 
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APPENDIX  A.  ESTIMATION  ALGORITHMS 


A.  NOTATION 

The  following  notation  is  a  generalization  of  that  used 
by  Robinson  [Ref.  2:p.  71]  as  it  pertains  to  the  aggregation 
scheme  used  in  this  paper.  The  indexing  system  has  been 
changed  to  reflect  the  use  of  officer  grades  to  identify  the 
cells . 

i  =  index  of  MOS  cells  in  the  aggregate, 
j  =  index  of  LOS  cells  in  the  aggregate, 
k  =  index  of  Grade  cells  in  the  aggregate, 
t  -  index  of  time  periods  in  the  aggregate. 
inv(i,j,k,t)  =  inventory  with  MOS  i,  LOS  j  and  Grade  k 
during  year  t. 

y(i,j,k,t)  »  number  of  attritions  in  cell  (i,j,k)  during 
year  t. 

n(i,j,k,t)  =  maximum  I y ( i , j , k , t ) , . 5 [inv ( i , j , k , t )  + 
inv (i,j,k,t  +  1 ) ] I 

D  =  matrix  identifying  cells  with  no  inventory  over  all 
estimation  years. 

d(i,j,k)  =  0,  if  cell  is  a  structural  zero. 
d(i,j,k)  =  1,  if  cell  is  not  a  structural  zero. 
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B.  FREEMAN-TUKEY  EXACT  INVERSION  ALGORITHM 

This  algorithm  provides  a  formula  for  the  inverse  of  the 
Freeman-Tukey  double  arcsine  transformation  to  express  the 
means  of  the  transformed  values  into  proportions  on  the 
original  scale. 

STEP  1.  Transform  the  data  using  the  Freeman-Tukey 
double  arcsine  transformation. 


x  (i  ,  j  ,  k ,  t )  =  0.5[n(i,j,k,t)+0.5]l/2  {sin1  [  L 

(n ( 1 , 3 , k , t ) +1 ) 


■1] 


+ 


sin*  1 


r2 (y (i , j , k, t) +1 
L (n (i , j , k , t ) +1 ) 


-1]  I 


^  \ 


STEP  2.  Calculate  the  average  cell  inventory  and 
transformed  values  over  all  T. 


x  (i  ,  j  ,k)  =  (1/T)Z!  x  ( i  ,  j , k)  ,  (A. 2) 

t 


inv  (i  ,  j  ,  k)  =  (1/T)Z:  inv(i,j,k).  (A. 3) 

t 


STEP  3.  Check  for  inversion  values  less  than  zero 
or  greater  than  one.  The  algorithm  is  terminated  at  this 
point  if  these  conditions  are  met. 

tm  =  x  (i  ,  j  ,  k)  /  ( inv  (i ,  j  ,  k)  +0 . 5) 1  /  2  +  n/2,  (A. 4) 

tmi  =  sin*  1  (1/  (inv(i ,  j  ,  k)  +  DM'*,  (A. 5) 


If 

tm  < 

tmi  , 

FTE  =  0.0, 

IF 

tm  > 

(n-tmi ) , 

FTE  =  1.0, 
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STEP  4.  Compute  the  sgn  function. 


If  cos(tm) 

>  0.0 

sgn  =*  1.0, 

If  cos { tm) 

<  0.0 

sgn  =  -1.0. 

STEP  5.  Computer  the  inverted  value  of  the 
transformed  variable 


FTE 


0 . 5 ( 1-sgn 1 1- [sin ( tm) 


sin(tm)  -  l/sinttm)-.  2  .  , 
inv ( i , j , k) 


(A.  6) 


C.  TRANSFORM  SCALE  VARIANCE  ALGORITHM 

This  algorithm  produces  the  variance  of  the  transformed 
data  within  the  unstable  variance  range  where  the  assumption 
of  normality  does  not  hold.  As  discussed  in  Chapter  II. 
this  range  of  instability  was  determined  to  be  between  1.001 
and  2.2.  The  method  for  computing  the  variance  within  this 
range  is  presented  here. 

STEP  1.  Compute  the  transformed  scale  values  to  be 
used  for  the  variance.  Since  x(i,j,k)  was  computed  in  the 
FTE  inversion  algorithm  of  this  appendix,  the  formulae  will 
not  be  repeated. 

x ( t )  =  21  [0.5  +  inv ( i , j , k , t )] 1 ' 2  ,  (A. 7) 

ijk 

z(i,j,k)  =  x(i,j,k)  +  (n/2)  .  (A.  3) 
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STEP  2.  Compute  the  variance  of  transformed  cell 
values.  The  formula  is  derived  in  Chapter  II  as  a  linear 
regression  equation.  The  range  of  the  variance  is  limited 
to  0.05  <  var(i,j,k)  5  1.0. 

var { i , j , k)  =  az(i,j,k)bl(z(i,j,k)  -  1 ) b  2  , 

1.001  <.  z  (i  ,  j  ,  k)  2  2.2  (A.  9) 

where  a  =  1.6835,  bl  =  -0.3934,  and  b2  =  0.3991. 

D.  EMPIRICAL  BAYES  SHRINKAGE  ALGORITHM 

This  algorithm  computes  the  empirical  Bayes  shrinkage 
values  and  estimator  of  attrition  rates.  Prior  to  entering 
this  algorithm,  the  value  for  a  has  been  set  to  0.0. 

STEP  1.  Initialization. 

ai  =  a .  ( A. 10 ) 

STEP  2.  Define  the  a  and  7*  parameters  in  the 
empirical  Bayes  context. 

a  =  ZI  l/(a  +  var  ( i  ,  j  , k)  )  ,  (A.  11) 

ijk 

where  var(i,j,k)  is  defined  by  formula  (A. 9). 

r  =  SI  a/  2  1/  (a  +  var  (  i  ,  j  ,  k )  )  .  (A. 12 

ijk  ijk 
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STEP  3.  Compute  the  empirical  Bayes  prior  variance. 


k-1-  23  a[x(i,  j  ,k)  -z]2 

a  =  a  -  _ ijk _ _  ,  (A. 13) 

21  [a  (x(i ,  j  ,  k)  -  z]  = 

ijk 


where  z  =  23  7*x(i,j,k)  and  k  =  23  d(i,j,k). 
ijk  ijk 


STEP  4.  Check  for  the  following  conditions  and 
branch  accordingly. 


If.  a  <  0.0  set  a  =  0.0  and  go  to  step  5. 

If,  |a  -  ai |  >  .0001  return  to  Step  1. 


STEP  5.  Compute  the  empirical  Bayes  estimator. 


::E  8  ( i  ,  j  ,  k] 


r - § -  1  v 

a+var (i , j , k) 

r  var(i,fj,k) 
a+var ( i , j , k) J 


k) 


+ 


(A. 14) 


STEP  6.  Compute  the  empirical  Bayes  shrinkage  set. 

shrEB(i,j,k)  =  var  ( i  ,  j  ,  k)  /  (  a  +■  var(i,j,k)).  (A. 15) 

E.  MODIFIED  EMPIRICAL  BAYES  SHRINKAGE  ALGORITHM 

This  algorithm  differs  only  slightly  from  that  of  the 
empirical  Bayes  shrinkage  algorithm  listed  above,  and  for 
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that  reason,  only  the  pertinent  steps  will  be  displayed 
here . 

Consider  a  particular  cell  and  its  transformed  value 
x(t)  for  t=l , . . . , T  in  the  estimation  set.  To  form  a 
weighted  time  average,  set 

xt(t)  =  x(t)/[0.5  +  inv(t)]1/2  (A. 16) 

This  leads  to  a  modification  of  step  1  in  the  following 
form. 

STEP  1.  Compute  the  weighted  transformed  scale 

average . 

xt  =  =  S  xt(t)  (A. 17) 

t 

Since  the  weights  will  modify  the  variance  of  the  weighted 
transformed  scale  values. 

var[xt(t)]  =  var (i , j ,k) / (0. 5  +  inv(t))  (A. 18) 

var[xt(t)]  =  (i)2  S  var[xt(t)]  =  vt/T  (A. 19) 

t 

The  extreme  right  of  formula  A. 19  serves  to  define  vt , 
the  approximation  for  the  weighted  variance  in  the 
validation  set  for  that  cell. 

The  values  of  xt  and  var(xt)  are  carried  throughout  this 
version  of  empirical  Bayes  estimation  to  form  the  values 
XTEB  for  the  empirical  Bayes  attrition  rate.  For  the 
purposes  of  inversion  to  the  original  scale,  the  usual  Basic 
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formula  is  modified  since  the  weighting  factor,  0.5  +  inv, 
is  already  incorporated  into  the  transform  value.  The  new 
Basic  inversion  for  this  version  is  simply  as  follows: 

p  =  1/2 [1  +  sin (XTEB ) ]  (A. 20) 
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APPENDIX  B.  MARINE  CORPS  PERSONNEL  INVENTORY 
ATTRITION  ANALYSIS  PROGRAM 
PROGRAM  MARCOR 


MARINE  CORPS  PERSONNEL  INVENTORY  ANALYSIS  PROGRAM  * 

* 

PURPOSE:  TO  PROVIDE  FUTURE  PERSONNEL  LOSS  ESTIMATION  OF  THE  * 
MARINE  CORPS  OFFICER  STRUCTURE  BY  SEVERAL  DIFFERENT  * 
ESTIMATION  SCHEMES.  * 

PROGRAMMERS:  LUIS  URIBE,  INDEPENDENT  CONTRACTOR  * 

CAPT  C  R  DICKINSON,  USMC  * 

DESCRIPTION:  WRITTEN  IN  FORTRAN  77  FOR  THE  IBM  3033  MAINFRAME  * 
COMPUTER  RESIDENT  AT  NPS,  THIS  PROGRAM  INCORPORATES  * 
METHODS  EXPLORED  BY  MAJOR  D. D.  TUCKER  ,  MAJOR  J.  R.  * 
ROBINSON  AND  CAPTAIN  DICKINSON.  ESTIMATION  SCHEMES  * 

INCLUDED  ARE:  * 

MAXIMUM  LIKELIHOOD  ESTIMATION  (MLE)  * 

TRANSFORMED  SCALE  CELL  AVERAGE  (TSCA)  * 

JAMES -STEIN  * 

LIMITED  TRANSLATION  JAMES-STEIN  * 

EMPIRICAL  BAYES  ESTIMATION  * 

RISK  ANALYSIS  WAS  PERFORMED  IN  BOTH  THE  ORIGINAL  AND  * 


TRANSFORMED  SCALES.  ON  THE  TRANSFORMED  SCALE,  * 
ACHIEVED  BY  USING  THE  FREEMAN -TUKEY  DOUBLE  ARCSINE  * 
TRANSFORMATION,  THE  RISK  WAS  DEFINED  BY  THE  AVERAGE  * 
SQUARED  DEVIATION  OF  ACTUAL  AND  ESTIMATED  VALUES.  * 
FOR  THE  ORIGINAL  SCALE,  THE  TRANSFORMED  VALUES  WERE  * 
FIRST  INVERTED  USING  TWO  TECHNIQUES.  AN  AD-HOC  * 
FORMULA  USED  BY  TUCKER  AND  ROBINSON  AND  AN  EXACT  * 


INVERSE  FORMULA  FROM  JOHN  J  MILLER.  TWO  MEASURES  OF  * 
EFFECTIVENESS,  CHI-SQUARE  GOODNESS  OF  FIT  AND  MEAN  * 
ABSOLUTE  DEVIATION  (MAD),  WERE  THEN  EVALUATED.  * 
RESULTS  ARE  PRINTED,  WITH  ALL  ARRAYS  COMPUTED,  AT  THE  * 


COMPLETION  OF  THE  PROGRAM.  * 

VARIABLE  DESCRIPTION:  * 

THE  FOLLOWING  LIST  OF  VARIABLES,  WITH  A  BRIEF  * 

DESCRIPTION,  FORM  THE  BASIS  FOR  THE  VARIABLE  NAMES  * 
FOR  THE  ENTIRE  PROGRAM.  * 

MLE  -  MAXIMUM  LIKELIHOOD  ESTIMATOR  * 

TSCA  -  TRANSFORMED  SCALE  CELL  AVERAGE  * 

J  -  JAMES-STEIN  ESTIMATOR  * 

J1  -  JAMES-STEIN  ESTIMATOR  MODIFIED  BY  * 

A  VARIATION  ON  SHRINKAGE  * 

JLT  -  LIMITED  TRANSLATION  JAMES-STEIN  * 

EB  -  EMPIRICAL  BAYES  * 

TEB  -  EMPIRICAL  BAYES  MODIFIED  BY  A  * 

VARIATION  ON  * 

TRANSFORMED  SPACE  VARIABLES  * 

THE  FOLLOWING  ARE  UNIQUE  VARIABLES  IN  THE  TRANSFORMED  * 
SPACE:  * 

X  -  VALUE  OF  THE  FREEMAN-TUKEY  * 


MAR00010 

MAR00020 

MAR00030 

MAR00040 

MAROOOSO 

MAR00060 

MAR00070 

MAR00080 

MAR00090 

MAR00100 

MAR00110 

MAR00120 

MAR00130 

MAR00140 

MAR00150 

MAR00160 

MAR0OI7O 

MAR00180 

MAR00190 

MAR00200 

MAR00210 

MAR00220 

MAR00230 

MAR00240 

MAR00250 

MAR00260 

MAR00270 

MAR00280 

MAR00290 

MAR0C300 

MAR00310 

MAR00320 

MAR00330 

MAR00340 

MAR00350 

MAR00360 

MAR00370 

MAR00380 

MAR00390 

MAR00400 

MAR00410 

MAR00420 

MAR00430 

MAR00440 

MAR00450 

MAR00460 

MAR00470 

MAR00480 

MAR00490 
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★ 

TRANSFORMATION 

* 

MAR00500 

* 

XB  -  (X-BAR)  AVERAGE  TRANSFORM  VALUE 

* 

MAR00510 

* 

OVER  ESTIMATION  YEARS 

* 

MARG0520 

* 

XBB  -  GRAND  MEAN  OF  TRANSFORMED  VALUES 

MAR0Q530 

h 

SSE  -  SUM  OF 

SQUARED  DIFFERENCES,  ERROR 

* 

*  MAR00540 

* 

(X  AND 

X-BAR) 

* 

MAR00550 

* 

SST  -  SUM  OF 

SQUARED  DIFFERENCES,  TOTAL 

* 

MAR00560 

* 

(X  AND 

XBB) 

* 

MAR00570 

# 

SSB  -  SUM  OF 

SQUARED  DIFFERENCES,  BETWEEN 

* 

MAR00580 

* 

( SST-SSE) 

* 

MAR00590 

* 

SHR  - 

JAMES-STEIN  SHRINKAGE  VALUES 

* 

MAR00600 

* 

SHR1  - 

MODIFIED  JAMES-STEIN  SHRINKAGE 

* 

MAR00610 

* 

VALUES 

* 

MAR00620 

* 

SHREB  - 

EMPIRICAL  BAYES  SHRINKAGE  VALUES 

* 

MAR00630 

* 

SHRTEB  - 

MODIFIED  EMPIRICAL  BAYES  SHRINKAGE 

* 

MAR00640 

* 

VALUES 

* 

MAR00650 

* 

THE  FOLLOWING  ARE  MODIFIED  BASIC  VARIABLES  IN  THE 

* 

MAR00660 

* 

TRANSFORMED  SPACE. 

THIS  IS  DONE  BY  PLACING  AN  'X'  FOR 

* 

MAR00670 

* 

A  PURE  TRANSFORMED 

VARIABLE,  OR  AN  'R'  FOR  A  RISK 

* 

MAR00680 

* 

VARIABLE  IN  FRONT  v,F  THE  BASIC  VARIABLE. 

* 

MAR00690 

* 

MLE 

RML  (MLE) 

* 

MAR00700 

* 

TSCA 

RTSCA 

★ 

MAR00710 

* 

XJ 

RJ 

* 

MAR00720 

* 

XJ1 

RJ1 

* 

MAR00730 

* 

XJLT 

RJLT 

* 

MAR00740 

* 

XEB 

RXEB 

* 

MAR00750 

* 

XTEB 

RTEB 

* 

MAR00760 

* 

ORIGINAL  SPACE  VARIABLES 

* 

MAR00770 

* 

THE  FOLLOWING  ARE  MODIFIED  BASIC  VARIABLES  IN  THE 

* 

MAR00780 

* 

ORIGINAL  SPACE.  SIMILAR  TO  WHAT  IS  DONE  IN  THE 

* 

MAR00790 

* 

TRANSFORMED  SPACE, 

A  'P'  SIGNIFIES  A  PURE  INVERTED 

* 

MAR00800 

* 

VARIABLE,  AND  AN  'RR'  REPRESENTS  A  RISK  VARIABLE.  AN 

* 

MAR00810 

* 

'M‘  AT  THE  END  OF  A  VARIABLE  NAME  IDENTIFIES  A  MILLER 

* 

MAR00820 

* 

INVERSE  VARIABLE. 

* 

MAR00830 

* 

PMLE 

RRML  (MLE) 

* 

MAR00840 

* 

PJ 

RRJ 

* 

MAR00850 

* 

PJM 

RRPJM 

* 

MAR00860 

* 

PJ1 

RRJ1 

* 

MAR00870 

* 

PJ1M 

RRJ1M 

* 

MAR00880 

* 

PJLT 

ERJLT 

* 

MAR00890 

* 

PEB 

RRXEB 

* 

MAR00900 

* 

PEBM 

RRXEBM 

★ 

MAR00910 

* 

PTEB 

RRRTEB 

* 

MAR00920 

* 

PTEBM 

RRTEBM 

* 

MAR00930 

* 

THE  FOLLOWING  LIST 

IS  OF  PREFIXES  USED  TO  DESCRIBE 

* 

MAR00940 

* 

THE  MEAN  ABSOLUTE  DEVIATION  VARIABLES: 

* 

MAR00950 

* 

MO  -  MEAN  DEVIATION  ’OVERAGE' 

* 

MAR00960 

* 

MU  -  MEAN  DEVIATION  'UNDERAGE* 

* 

MAR00970 

* 

MD  -  MEAN  'TOTAL  DEVIATION 

* 

MAR00980 

* 

KVO  -  NUMBER 

OF  'OVERAGE'  ESTIMATES  IN  THE 

* 

MAR00990 

* 

AGGREGATE 

* 

MAROIOOO 

* 

KVU  -  NUMBER 

OF  'UNDERAGE*  ESTIMATES  IN  THE 

* 

MAR01010 

•k 

AGGREGATE 

* 

MAR01020 

* 

KVD  -  NUMBER 

OF  'TOTAL*  ESTIMATION  DEVIATIONS 

* 

MAR01030 

* 

IN  THE 

AGGREGATE 

* 

MAR01040 

* 

THESE  PREFIXES  ARE 

APPENDED  TO  THE  BASIC  VARIABLES  TO 

* 

MAR01050 
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* 

* 

* 

★ 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


FORM  THE  MAD  VARIABLES  AS  FOLLOWS  * 

MOTSCA  MUTSCA  MDTSCA  KVOTSA  KVUTSA  KVDTSA  (TSCA)  * 

MOTSAM  MUTSAM  MDTSAM  KVOTSM  KVUTSM  KVDTSM  (TSCAM)  * 

MOPMLE  MUPMLE  MDPMLE  KVOPML  KVUPML  KVDPML  (PMLE)  * 

MOPJ  MUPJ  MDPJ  KVOPJ  KVUPJ  KVDPJ  (PJ)  * 

MOPJM  MUPJM  MDPJM  KVOPJM  KVUPJM  KVDPJM  (PJM)  * 

MOPJ1  MUPJ1  MDPJ1  KVOPJ 1  KVUPJ1  KVDPJ1  (PJ1)  * 

MOPJ1M  MUPJ1M  MDPJ1M  KVOP1M  KVUP1M  KVDP1M  (PJ1M)  * 

MOPEB  MUPEB  MDPEB  KVOPEB  KVUPEB  KVDPEB  (PEB)  * 

MOPEBM  MUPEBM  MDPEBM  KVOPBM  KVUPBM  KVDPBM  (PEBM)  * 

MOPTEB  MUPTEB  MDPTEB  KVOPTB  KVUPTB  KVDPTB  (PTEB)  * 

MOPTBM  MUPTBM  MDPTBM  KVOBTM  KVUBTM  KVDBTM  (PTEBM)  * 


INPUT/OUTPUT:  * 

INPUT  IS  TAKEN  FROM  A  FILE  CREATED  FROM  A  TAPE  * 

PROVIDED  BY  NPRDC  * 

OUTPUT  IS  TO  THE  PRINTER  * 


**  FIXED  PARAMETERS 

PARAMETER  (MXYR=7,  MXLOS=31,  MXGRD=10,  MXMOS=15) 

**  INPUT  PARAMETERS 

INTEGER  ST1,  ST2,  LYR 
INTEGER  SLOS1,  SLOS2 
INTEGER  SMOS(MXMOS),  NMOS 
INTEGER  SGRD(MXGRD) ,  NGRD 
**  VARIABLES  IN  ORIGINAL  SPACE 

INTEGER  T,  FLAG,  IVYR(MXYR) 

INTEGER*2  D(MXMOS,  MXLOS,  MXGRD) 

INTEGER  KV(MXYR) 

REAL  AVCINV , AINV , ATT , ATT1 

REAL  CINV(MXMOS, MXLOS, MXGRD, MXYR) ,Y(MXMOS , MXLOS, MXGRD, MXYR) 
REAL  AC INV( MXMOS, MXLOS, MXGRD) 

REAL  PJ(MXMOS, MXLOS, MXGRD),  PJMCMXMOS, MXLOS, MXGRD) 

REAL  PMLE(MXMOS, MXLOS, MXGRD),  PJLT( MXMOS, MXLOS, MXGRD) 

REAL  PJl(MXMOS, MXLOS, MXGRD),  PJ1M(MXM0S, MXLOS, MXGRD) 

REAL  PEB (MXMOS, MXLOS, MXGRD ) ,  PEBM(MXMOS, MXLOS, MXGRD) 

REAL  PTEB ( MXMOS , MXLOS , MXGRD ) ,  PTEBMC MXMOS , MXLOS , MXGRD ) 

REAL  RRTSCA(MXYR) .RRJ(MXYR) .RRJLT(MXYR) .RRML(MXYR) 

REAL  RRTSAM(MXYR) .RRJM(MXYR) .RRJl(MXYR) ,RRJ1M(MXYR) 

REAL  RRPEB(MXYR) ,RRPEBM(MXYR) 

REAL  RRPTEBC MXYR ),RRPTBM( MXYR) 

**  VARIABLES  IN  TRANSFORMED  SPACE 

REAL  X( MXMOS, MXLOS, MXGRD, MXYR) 

REAL  XT( MXMOS, MXLOS, MXGRD, MXYR) 

REAL  XBCMXMOS, MXLOS, MXGRD),  TSCA( MXMOS, MXLOS, MXGRD) 

REAL  XTB(MXMOS, MXLOS, MXGRD) 

REAL  X J( MXMOS , MXLOS , MXGRD ) ,  X JLT( MXMOS , MXLOS , MXGRD ) 

REAL  XEB( MXMOS, MXLOS, MXGRD),  XJl(MXMOS .MXLOS , MXGRD) 

REAL  XTEB( MXMOS, MXLOS, MXGRD) 

REAL  XMLE( MXMOS, MXLOS, MXGRD) 

REAL  SHREB(MXMOS, MXLOS, MXGRD) 

REAL  SHRTEB( MXMOS, MXLOS, MXGRD) 

REAL  RML(MXYR) ,  RSL(MXYR) .RXEB(MXYR) 

REAL  RXTEB(MXYR) 

REAL  RTSCA( MXYR ) , RTJ( MXYR ) , RTJLTC MXYR ) 

REAL  RJ( MXYR), RJLT( MXYR),  RJl(MXYR) 

REAL  TSCAM( MXMOS, MXLOS, MXGRD) 


MAR01060 
MAR01070 
MAR01080 
MAR01090 
MAROllOO 
MAROlllO 
MAR01 120 
MAR01130 
MAR01140 
MAR01150 
MAR01160 
MAR01170 
MAR01180 
MAR01190 
MAR01200 
MAR01210 
MAR01220 
MAR01230 
MAR01240 
MAR01250 
MAR01260 
MAR01270 
MAR01280 
MAR01290 
MAR01300 
MAR01310 
MAR01320 
MAR01330 
MAR01340 
MAR01350 
MAR01360 
MAR01370 
MAR01380 
MAR01390 
MAR01400 
MAR01410 
MAR01420 
MAR01430 
MAR01440 
MAR01450 
MAR01460 
MAR01470 
MAR01480 
MAR01490 
MAR01500 
MAR01510 
MAR01520 
MAR01530 
MAR01540 
MAR01550 
MAR01560 
MAR01570 
MAR01580 
MAR01590 
MAR01600 
MAR01610 
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REAL  V( MXMOS , MXLOS , MXGRD ) 

REAL  VT( MXMOS, MXLOS, MXGRD) 

REAL* 8  SHR,  SHRX,  SKR1,  SSB,  SST,  SSE,  XBB 
**  VARIABLES  FOR  FIGURE  OF  MERIT 

INTEGER  KVOTSA(MXYR) .KVOTSM(MXYR) .KVOPJ(MXYR) .KVOPJl(MXYR) 

INTEGER  KVOPJM(MXYR) .KVOPEB(MXYR) ,KV0P1M(MXYR) .KVOPBM(MXYR) 
INTEGER  KVOPTB(MXYR) ,KVOBTM(MXYR) 

INTEGER  KVOPML(MXYR) 

INTEGER  KVUTSA(MXYR) .KVUTSM(MXYR) .KVUPJ(MXYR) .KVUPJl(MXYR) 

INTEGER  KVUPJM(MXYR) .KVUPEB(MXYR) .KVUPIM(MXYR) .KVUPBM(MXYR) 
INTEGER  KVUPTB(MXYR) .KVUBTM(MXYR) 

INTEGER  KVUPML(MXYR) 

INTEGER  KVDTSA(MXYR) .KVDTSM(MXYR) .KVDPJ(MXYR) .KVDPJl(MXYR) 

INTEGER  KVDPJM(MXYR) ,KVDPEB(MXYR) .KVDPIM(MXYR) ,KVDPBM(MXYR) 
INTEGER  KVDPTB(MXYR) .KVDBTM(MXYR) 

INTEGER  KVDPML(MXYR) 

REAL  MUTSCA(MXYR) .MUTSAM(MXYR) .MUPJ(MXYR) .MUPJl(MXYR) .MUPJM(MXYR) 
REAL  MUPJIM(MXYR) .MUPEBM(MXYR) .MUPEB(MXYR) .MUPMLE(MXYR) 

REAL  MUPTBM(MXYR) .MUPTEB(MXYR) 

REAL  MOTSCA(MXYR) ,MOTSAM(MXYR) .MOPJ(MXYR) .hOPJl(MXYR) .MOPJM(MXYR) 
REAL  MOPJIM(MXYR) .MOPEBM(MXYR) .MOPEB(MXYR) ,MOPMLE(MXYR) 

REAL  MOPTBM(MXYR) ,MOPTEB(MXYR) 

REAL  MDTSCA(MXYR) .MDTSAM(MXYR) .MDPJ(MXYR) .MDPJl(MXYR) .MDPJM(MXYR) 
REAL  MDPJIM(MXYR) .MDPEBM(MXYR) ,MDPEB(MXYR) .MDPMLE(MXYR) 

REAL  MDPTBM(MXYR) .MDPTEB(MXYR) 

REAL  MILLER 

**  INPUT  DATA  AREAS 

INTEGER  MOS,  LOS,  GRADE,  YR,  INV 
INTEGER  MOS 1, LOS 1, GRADE 1,YR1,INV1 
**  INPUT  INITIALIZATION 

DATA  KV/MXYR*0/ 

DATA  MOPMLE/MXYR*0/,MUPMLE/MXYR*0/,MDPMLE/MXYR*0/ 

DATA  KVOPML/MXYR*0/ , KVUPML/MXYR*0/ ,  KVDPML/MXYR*0/ 

DATA  MOTSCA/MXYR*0/ ,MUTSCA/MXYR*0/ ,MDTSCA/MXYR*0/ 

DATA  KVOTSA/MXYR*0/ ,KVUTSA/MXYR*0/ ,KVDTSA/MXYR*0/ 

DATA  MOTSAM/MXYR*0/,MUTSAM/MXYR*0/ ,MDTSAM/MXYR*0/ 

DATA  KVOTSM/MXYR*O/,KVUTSM/MXYR*0/ , KVDTSM/MXYR*0/ 

DATA  M0PJ/MXYR*0/ ,MUPJ/MXYR*0/ ,MDPJ/MXYR*0/ 

DATA  KVOPJ/MXYR*0/ ,KVUPJ/MXYR*0/ ,KVDPJ/MXYR*0/ 

DATA  M0PJM/MXYR*0/ ,MUPJM/MXYR*0/ ,MDPJM/MXYR*0/ 

DATA  KVOPJM/MXYR*0/ ,KVUPJM/MXYR*0/ ,KVDPJM/MXYR*0/ 

DATA  MOPJ1/MXYR*0/ ,MUPJ1/MXYR*0/ ,MDPJ1/MXYR*0/ 

DATA  KVOPJ1/MXYR*0/,KVUPJ1/MXYR*0/,KVDPJI/MXYR*0/ 

DATA  MOPJ1M/MXYR*0/,MUPJ1M/MXYR*0/ ,MDPJ1M/MXYR*0/ 

DATA  KVOP1M/MXYR*0/ ,KVUP1M/MXYR*0/ ,KVDP1M/MXYR*0/ 

DATA  MOPEB/MXYR*0/ ,MUPEB/MXYR*0/ ,MDPEB/MXYR*0/ 

DATA  KVOPEB/MXYR*0/ ,KVUPEB/MXYR*0/ ,KVDPEB/MXYR*0/ 

DATA  MOPEBM/MXYR*0/ ,MUPEBM/MXYR*0/ , MDPEBM/MXYR*0/ 

DATA  KV0PBM/MXYR*0/,KVUPBM/MXYR*0/ ,KVDPBM/MXYR*0/ 

DATA  MOPTEB/MXYR*0/,MUPTEB/MXYR*0/ ,MDPTEB/MXYR*0/ 

DATA  KVQPTB/MXYR*O/,KVUPTB/MXYR*0/,KVDPTB/MXYR*0/ 

DATA  MOPTBM/MXYR*0/ ,MUPTBM/MXYR*0/ , MDPTBM/MXYR*0/ 

DATA  KVOBTM/MXYR*0/,KVUBTM/MXYR*0/,KVDBTM/MXYR*0/ 

Af*AAA>*A*A»*A*A»*A»*«r»**A  A  A  AAA  A  A*  A  A»*»**A*A*AA  A  A  AAA  A  A  A  AAAAAAAAA******* 


MAR01620 
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MAR01650 
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MAR01700 
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MAR01740 
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MAR01760 

MAR01770 

MAR01780 

MAR01790 

MAR01800 

MAR01810 

MAR01820 

MAR01830 

MAR01840 

MAR01S50 

MAR01860 

MAR01870 

MAR01880 

MAR01890 

MAR01900 

MAR01910 

MAR01920 

MAR01930 

MAR01940 

MAR01950 

MAR01960 

MAR01970 

MAR01980 

MAR01990 

MAR02000 

MAR02010 

MAR02020 

MAR02030 
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MAR02150 

MAR02160 
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**  READ  PARAMETER  DATA  TO  USE  IN  THIS  RUN.  ** 

a********************************************************1************* 

WRITE(5.100) 

100  FORMATC '  ENTER  1ST  YEAR  AND  LAST  YEAR  TO  USE  FOR  ESTIMATION') 
READ(5 ,*)  ST1,ST2 

WRITE( 6 ,*)  'ESTIMATION  YEARS:  '.ST1.ST2 

* 

WRITEC5.101) 

101  FORMATC  ENTER  LAST  YEAR  PRESENT  IN  THE  DATA  BASE') 

R£AD(5 ,*)  LYR 

WRITE(6,*)  'LAST  YEAR  AVAILABLE  FOR  VALIDATION:  ' ,LYR 

* 

WRITE(5. 102) 

102  FORMATC1  ENTER  NO.  OF  MOS  FOLLOWED  BY  ARRAY  OF  SAME  LENGTH') 
READ(5 ,*)  NMOS,  (SMOS(I),  I=1,NM0S) 

WRITEC6,*)  'MOS  SELECTED:',  (SMOS(I),  I=1,NM0S) 

* 

WRITEC5.103) 

103  FORMATC f  ENTER  1ST  AND  LAST  LOS  VALUE  TO  USE’) 

READC5,*)  SL0S1 ,SL0S2 

WRITEC6,*)  'LOS  RANGE:  ’ .SL0S1.SL0S2 

* 

WRITE(5,104) 

104  FORMATC  ENTER  NO.  OF  GRADE  FOLLOWED  BY  ARRAY  OF  SAME  LENGTH’) 
READC5,*)  NGRD,  C SGRDC I ) ,  I=1,NGRD) 

WRITEC6,*)  'GRADES  SELECTED',  (SGRDCI),  1=1, NGRD) 

* 

WRITEC5.105) 

105  FORMATC  ENTER  DEE  FACTOR') 

READC5,*)  DEE 

WRITEC6,*)  'DEE  FACTOR: ', DEE 

* 

**  COMPUTE  START  AND  END  YEARS  FOR  VALIDATION 

NYR=LYR-ST1+1 
NVYR=LYR-ST2 
NL0S=SL0S2-SL0S1+1 
T=ST2-ST1+1 
DO  41  1=1 ,NVYR 
IVYRC I)=ST2+I 
41  CONTINUE 

*  *  *  *  *  A  A***  MriHrt  ****  A  ****  *****  A  *  A  *  **  *****  ************ 


**  INITIALIZE  D  ARRAY  TO  ALLOW  FINDING  MISSING  VALUES.  ** 

k*kk******A***k*A*/t**k*h*+lr****A****+*/,****+*.;\  A  A*  **************  ******* 

DO  1  I=1,MXM0S 
DO  2  J=1 ,MXLOS 
DO  3  K=1,MXGRD 
DC  I , J,K)=-9999 
3  CONTINUE 

2  CONTINUE 
1  CONTINUE 

irk***  A  it****  A  ****************** *********  *********  ********************** 

**  ■  TAPE  PROCESSING  TO  INPUT  INVENTORY  AND  ATTRITION  VALUES,  ** 

**  DEVELOPE  D  MATRIX  ACCORDING  TO  EXTRACTION  CRITERIA.  ** 

***********  A  ill  **A**A**ir**A**A***A  ******  A  ***  *********************  ******** 

FLAG=0 

IEOF=0 
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MAR02180 

MAR02190 
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MAR02330 

MAR02340 

MAR02350 

MAR02360 

MAR02370 

MAR02380 

MAR02390 

MAR02400 

MAR02410 

MAR02420 

MAR02430 

MAR02440 

MAR02450 

MAR02460 

MAR02470 

MAR02480 

MAR02490 

MAR02500 

MAR02510 

MAR02520 

MAR02530 

MAR02540 

MAR02550 

MAR02560 

MAR02570 

MAR02580 

MAR02590 

MAR02600 
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MAR02630 

MAR02640 

MAR02650 

MAR02660 

MAR02670 

MAR02680 
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* 


* 

11 


** 


** 

* 

12 

** 


8 

** 


CALL  R£ADER( MOS 1 , LOS 1 , GRADE 1 , YR 1 , I NV 1 , ATTl ,  NMOS.SMOS, 

*  NGRD.SGRD,  SLOS1.SLOS2,  ST1 , ST2 ,  IM1,IL1,IG1,IT1,IE0F) 

IF(IEOF  .NE.  0)  THEN 

WRITE(6,*)  '**★*  No  DATA  MEETS  SELECTIONS  REQS* 

STOP 

END  IF 

IF( IEOF  .NE.  0)  GO  TO  9 
MOS-MOS1 
LOS-LOS1 
GRADE-GRADE 1 
YR-YR1 
INV-INV1 
ATT-ATT1 
IM-IM1 
IL-IL1 
IG-IG1 
IT-IT1 

IF(YR. GT.  ST1)  FLAG  =  1 
CHECK  CASE  WHERE  DATA  BEGINS  PAST  1ST  YR 
IF(YR  .GT.  ST1)  CINV( IM , IL , IG , IT- 1 )=.  5*FLOAT( INV) 

CALL  READER ( MOS 1 , LOS  1 , GRADE 1 , YR 1 , I NV 1 , ATT 1 ,  NMOS.SMOS, 

*  NGRD.SGRD,  SLOS1.SLOS2,  ST1.ST2,  IM1,IL1,IG1,IT1,IE0F) 
ONLY  FOR  CHECKING  PURPOSES.  MARK  D  TO  INDICATE  SOME  YR  PRESENT 

D(IM,IL,IG)=0 


MOS1.EQ.MOS  .AND.  LOS1.EQ.LOS 
GRADE1.EQ.  GRADE))  GO  TO  8 
1 


IF(.  NOT.  (IEOF.  EQ.O  .AND. 

.AND. 

CENTRAL  INV.  FOR  YR1  - 
IF( YR1  .GT.  YR+1)  THEN 

CINVC IM1 , IL1 , IG1 , IT1- 1)=.  5*FLOAT( 1NV1) 
AINV=.  5*FLOAT( INV) 

I = YR+1 
FLAG=1 


MAR02730 

MAR02740 

MAR027S0 

MAR02760 

MAR02770 

MAR02780 

MAR02790 

MAR02800 

MAR02810 

MAR02820 

MAR02830 

MAR02840 

MAR028S0 

MAR02860 ’ 

MAR02870 

MAR02880 

MAR02890 

MAR02900 

MAR02910 

MAR02920 

MAR02930 

MAR02940 

MAR02950 

MAR02960 

MAR02970 

MAR02980 

MAR02990 

MAR03000 

MAR03010 

MAR03020 

MAR03030 

MAR03040 

MAR03050 

MAR03060 

MAR03070 


ELSE 

AINV-.  5*FLOAT( INV+INVl) 

END  IF 

CINVC IM, IL, IG, IT)=AMAX1( AINV , ATT) 

Y( IM, IL,IG, IT)=ATT 

YR-YR1 

INV-INV1 

ATT-ATT1 

IM-IM1 

IL-IL1 

IG-IG1 

IT-IT1 

CALL  READER(M0S1,L0S1,GRADE1,YR1,INV1,ATT1,  NMOS.SMOS, 

*  NGRD.SGRD,  SLOS1.SLOS2,  ST1.ST2,  IM1,IL1,IG1,IT1,IE0F) 

GO  TO  12 
CONTINUE 

WHEN  YEARS  MISSING  AT  THE  END 
IF(YR  .LT.  LYR)  THEN 
FLAG-1  •- 
AINV-. 5*FLOAT( INV) 

ELSE 


MAR03080 

MAR03090 

MAR03100 

MAR03110 

MAR03120 

MAR03130 

MAR03140 

MAR03150 

MAR03160 

MAR03170 

MAR03180 

MAR03190 

MAR03200 

MAR03210 

MAR03220 

MAR03230 

MAR03240 

MAR03250 

MAR03260 

MAR03270 

MAR03280 
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AINV=FLOAT( INV) 

END  IF 

CINV( IM, IL, IG, IT)*AMAX1(AINV,ATT) 
Y( IM,IL,IG,IT)=ATT 


CC=CINV( IM , IL , IG , IT) 

TEMP=- 1. +2. *ATT/( 1. +CC) 

TEMP1«-1.  +  2.  *(  1.  +ATT)/(  1.  +CC) 

COABS(TEMP) 

CC1=ABS( TEMPI) 

IF(CC  .GT.  1.  .OR.  CC1  .  GT.  1.)  VRITE(6,*)  '***  TEMP,  TEMPI*’ , 
*  TEMP, TEMP 1,'***M0S, LOS, GR,YR=* ,SM0S(IM) ,IL,SGRD(IG) ,IT 

* 

GO  TO  11 
9  CONTINUE 


** 

** 

****** 


CALCULATE  AVERAGE  CENTRAL  INVENTORY  AND  MAXIMUM  ** 

LIKELIHOOD  ESTIMATE.  ** 

******************************************* ********************* 


DO  10  IM=1 ,NMOS 
DO  20  IL=l,NLOS 
DO  30  IG=1,NGRD 
SUMY=0 

ACINV(IM,IL,IG)=0 
DO  35  IT=1,T 
SUMY=SUMY+Y( IM , IL , IG , IT) 

ACINV( IM, IL, IG)=ACINV( IM, IL, IG)+CINV( IM, IL, IG, IT) 

35  CONTINUE 

IF(ACINV( IM, IL,IG)  .  NE.  0)  THEN 

PMLE( IM , IL, IG)=SUMY/ACINV( IM, IL, IG) 

ELSE 

PMLE( IM,IL,1G)=Q 

END  IF 

AC INV ( IM, IL, IG)=ACINV( IM, IL, IG)/T 
DO  40  IT=1 ,NYR 

**  IF  ANY  CENTRAL  INVENT0RY>0  (OVER  EST.  YRS  T)  THEN  D=1 

IF(CINV(IM,IL,IG,IT).  GT.  0  .AND.  IT.  LE.  T)  D( IM, IL, IG)=1 
**  IF  ANY  CENTRAL  INVENT0RY>0  (OVER  VAL.  YRS  >  T)  THEN  COUNT 

IF(CINV(IM,IL,IG,IT).  GT.  0  .AND.  IT.  GT.  T)  KV( IT-T)=KV( IT-T)+1 
**  TRANSFORMATION  OF  CINV  USING  THE  FREEMAN-TUKEY 

**  DOUBLE  ARCS IN  FORMULA. 

CALL  FTT(CINV( IM, IL, IG, IT) ,Y( IM, IL,IG,IT) ,X( IM, IL, IG, IT) ) 
XT(IM,IL,IG,IT)=X(IM,IL,IG,IT)  /  SQRT(0. 5+CINV(IM,IL,IG,IT)) 
40  CONTINUE 

30  CONTINUE 
20  CONTINUE 
10  CONTINUE 

ft  A' *****  A  A  A  A  ArA  A  ill  A  »*■*■**  *■*;■*■**■*  A  ******* '*■*-*  *'*  A  A* r'A'"A"A  A  A  A  A  **** **************** 


**  CHECK  FOR  MISSING  DATA,  COMPUTE  THE  TRANSFORMED  ** 

**  CELL  AVERAGES  AND  GRAND  MEAN.  ** 

XBB-0 

KK=0 

DO  50  IM=l,NMOS  f* 

DO  60  IL=l,NLOS 
DO  70  IG*1,NGRD 


MAR03290 

MAR03300 

MAR03310 

MAR03320 

MAR03330 

MAR03340 

MAR03350 

MAR03360 

MAR03370 

MAR03380 

MAR03390 

MAR03400 

MAR03410 

MAR03420 

MAR03430 

MAR03440 

MAR03450 

MAR03460 

MAR03470 

MAR03480 

MAR03490 

MAR03500 

MAR03510 

MAR03520 

MAR03530 

MAR03540 

MAR03550 

MAR03560 

MAR03570 

MAR03580 

MAR03590 

MAR03600 

MAR03610 

MAR03620 

MAR03630 

MAR03640 

MAR03650 

MAR03660 

MAR03670 

MAR03680 

MAR03690 

MAR03700 

MAR03710 

MAR03720 

MAR03730 
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MAR03750 

MAR03760 

MAR03770 

MAR03780 

MAR03790 

MAR03800 

MAR03810 
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MAR03840 
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IF(D(IM,IL,IG)  .EQ.  -9999)  THEN 

MAR03850 

REPORT  MISSING  COMBINATION  IM.IL.IG  &  CLEAR  D  ENTRY 

MAR03860 

D(IM,IL,IG)=0 

MAR03870 

ILD=SL0S1+IL-1 

MAR03880 

FLAG=1 

MAR03890 

END  IF 

MAR03900 

KK*KK+D( IM.IL.IG) 

MAR03910 

XB( IM.IL.IG )=0 

MAR03920 

XTB(IM,IL,IG)=0 

MAR03930 

DO  80  IT=1,T 

MAR03940 

XB(IM,IL,IG)=XB(IM,IL,IG)+X( IM.IL.IG, IT) 

MAR03950 

XTB(IM,IL,IG)*XTB(IM,IL,IG)+XT( IM.IL.IG, IT) 

MAR03960 

80 

CONTINUE 

MAR03970 

XB(IM,IL,IG)=XB(IM,IL,IG)/T 

MAR03980 

XTB( IM, IL, IG)=XTB( IM.IL, IG)/T 

MAR03990 

** 

EMPIRICAL  BAYES  PREPARATION  ON  THE  TRANSFORMED  SCALE. 

MAR04000 

XEB( IM, IL, IG)“XB( IM, IL, IG) 

MAR04010 

XTEB(IM,IL,IG)=XTB( IM.IL.IG) 

MAR04020 

* 

MAR04030 

XBB=XBB+D( IM, IL, IG)*XB( IM , IL, IG) 

MAR0  040 

70 

CONTINUE 

MAR04050 

60 

CONTINUE 

MAR04060 

50 

CONTINUE 

MAR04070 

XBB-XBB/KK 

MAR04080 

MAR04090 

** 

BAYES  VARIANCE  ITERATION  SECTION.  ** 

MAR04100 

MAR04110 

IC=0 

MAR04120 

** 

COMPUTE  TRANSFORMED  SCALE  CELL  VARIANCE. 

MAR04130 

499 

DO  500  IM=1,NM0S 

MAR04140 

DO  510  IL=1,NL0S 

MAR04150 

DO  520  IG=1,NGRD 

MAR04160 

IF  (D(IM.IL.IG).  NE.  1)  GOTO  520 

MAR04170 

AVSQR  ==0.0 

MAR04180 

DO  522  IT=1,T 

MAR04190 

522 

AVSQR  ■  AVSQR  +  SQRT(0.5  +  CINV(IM,IL,IG,IT)) 

MAR04200 

Z  *  XEB(IM.IL.IG)  +  (AVSQR/T)  *  1.5708 

MAR04210 

V( IM.IL.IG)  =  VAR(Z) 

MAR04220 

520 

CONTINUE 

MAR04230 

510 

CONTINUE 

MAR04240 

500 

CONTINUE 

MAR04250 

* 

MAR04260 

A-Q 

MAR04270 

530 

AI*A 

MAR04280 

SALPH=0 

MAR 04290 

DO  531  IM=1,NM0S 

MAR04300 

DO  532  IL=1 ,NLOS 

MAR04310 

DO  533  IG=1,NGRD 

MAR04320 

533 

IF  ( D( IM , IL , IG) .  EQ.  1 )  SALPH  =  SALPH  +  1. /(A  +  V(IM,IL,IG)) 

MAR04330 

532 

CONTINUE 

MAR04340 

531 

CONTINUE 

MAR04350 

it 

4 

MAR04360 

ZB=0 

MAR04370 

irk 

COMPUTE  BAYES ’ALPHA  AND  GAMMA  PARAMETERS 

MAR04380 

DO  537  IM*1,NM0S 

MAR04390 

DO  538  IL=1,NL0S 

MAR04400 
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DO  539  IG=1,NGRD 

IF  (D(IM,IL,IG).EQ.  1)  THEN 
ALPH  =  l./(A  +  V(IM,IL,IG)) 

GAM  =  ALPH  /  SALPH 

ZB  ■  ZB  +  GAM*XB(IM,IL,IG) 

END  IF 

539  CONTINUE 
538  CONTINUE 
537  CONTINUE 

F=0 

G=0 

DO  540  IM=1 ,NMOS 
DO  541  IL=1 ,NLOS 
DO  542  IG=1,NGRD 

IF  (D(IM,IL,IG).  EQ.  1)  THEN 
ALPH  =  l./(A  +  V(IM,IL,IG>) 

F  =  F  +  (ALPH  *  (XBCIM.IL, IG)  -  ZB)**2) 

G  ■  G  +  (ALPH  *  (XB( IM, IL, IG) -ZB) )**2 

END  IF 

542  CONTINUE 

541  CONTINUE 

540  CONTINUE 

**  COMPUTE  THE  PRIOR  VARIANCE  (A) 

A  =  A  -  (KK-1-F)/G 

* 

IF  (A.LE.O.  )  THEN 
A=0 

GOTO  550 
ENDIF 

IF  (ABS(A-AI).GT.  .0001)  GOTO  530 
**  COMPUTE  THE  EMPIRICAL  BAYES  ATTRITION  RATES 

550  DO  570  IM=1,NM0S 

DO  580  IL=1,NL0S 
DO  590  IG»1,NGRD 

IF  (D(IM.IL.IG).EQ. 1)  THEN 
DEN  =  A  +  V( IM, IL, IG) 

XEB( IM, IL, IG)  =  (A/DEN)*XB( IM, IL, IG)+( (DEN-A)/DEN  )*ZB 
ENDIF 

590  CONTINUE 

580  CONTINUE 

570  CONTINUE 
IC=IC+1 

IF  (A.EQ.  0)  GOTO  599 
IF  (IC.  LT.  10)  GOTO  499 

★ 

599  CONTINUE 
* 

**  COMPUTE  THE  EMPIRICAL  BAYES  SHRINKAGE  RATES 

DO  600  IM=1,NM0S 
DO  601  IL=1 ,NLOS 
DO  602  IG=1,NGRD 
SHREBC IM, IL, IG)  =  0 
IF  (D(IM,IL,IG).  EQ.  1)  THEN 

SHREB( IM,IL,IG)  =  V( IM, IL, IG)/(A  +  V(IM,IL,IG)) 

ENDIF 

602  CONTINUE 
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CONTINUE 

CONTINUE 


ALTERNATE  BAYES  VARIANCE  ITERATION  SECTION.  ** 

**********  ft****  *  *■*■*  **  ******************************************** 

COMPUTE  TRANSFORMED  SCALE  CELL  VARIANCE. 

DO  700  IM=1 ,NMOS 
DO  710  IL=1 ,NLOS 
DO  720  IG=1,NGRD 
VT( IM, IL, IG)  =0.0 
IF  (D(IM,IL,IG).NE.  1)  GOTO  720 
DO  725  IT-1, T 

ZT  =  X(IM,IL,IG,m  +  1-  5708  *  SQRT(0.5  + 

*  CINV( IM, IL, IG , IT) ) 

VT( IM, IL, IG)  =  VT( IM,IL,IG)  +  VAR(ZT)/(0.5  + 

*  CINV( IM,IL,IG,IT)) 

VT(IM,IL,IG)  =  VT( IM, IL, IG)/T**2 

CONTINUE 

CONTINUE 

CONTINUE 

AT=0 

ATI=AT 

SALPHT=0 

DO  731  IM=1 ,NMOS 
DO  732  IL=1 ,NLOS 
DO  733  IG=1 ,NGRD 

IF  (D(IM,IL,IG).EQ. 1)  SALPHT=  SALPHT+  l./(AT+  VT(IM,IL,IG)) 
CONTINUE 
CONTINUE 
ZTB=0 

COMPUTE  BAYES  ALPHA  AND  GAMMA  PARAMETERS 
DO  737  IM=1 ,NMOS 
DO  738  IL=1,NL0S 
DO  739  IG=1 ,NGRD 

IF  (D(IM,IL,IG).  EQ.  1)  THEN 
ALPHT=  1.  /(AT+  VT( IM,IL,IG)) 

GAMT  =  ALPHT  /  SALPHT 

ZTB  =  ZTB+  GAMT*XTB(IM,IL,IG) 

ENDIF 

CONTINUE 

CONTINUE 

CONTINUE 

PZTB=0.  5*( 1.  +SIN( ZTB) ) 

FT=0 

GT=0 

DO  740  IM=1 ,NMOS 
DO  741  IL=1,NL0S 
DO  742  IG=1 ,NGRD 

IF  (D( IM, IL, IG). EQ. 1)  THEN 
ALPHT=  1. /(AT  +  VT(IM.IL.IG)) 

FT  =  FT  +  (ALPHT  *  (XTB( IM.IL.IG)  -  ZTB)**2) 

GT  =  GT  +  (ALPHT  *  (XTB( IM, IL, IG)  -  ZTB))**2 
ENDIF  '* 

CONTINUE 

CONTINUE 
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MAR04980 

MAR04990 

MAR05000 

MAR05010 

MAR05020 

MAR05030 

MAR05040 

MAR05050 

MAR05060 

MAR05070 

MAR05080 

MAR05090 
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MAR05110 

MAR05120 
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MAR05140 

MAR05150 
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MAR05180 
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MAR05280 

MAR05290 
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MAR05330 

MAR05340 

MAR05350 

MAR05360 

MAR05370 

MAR05380 

MAR05390 

MAR05400 

MAR05410 

MAR05420 

MAR05430 

MAR05440 

MAR05450 
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MAR05470 

MAR05480 

MAR05490 

MAR05500 
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740 

CONTINUE 

MAR05530 

** 

COMPUTE  THE  PRIOR  VARIANCE  (A) 

MAR05540 

AT=  AT-  (KK-1-FT)/GT 

MAR05550 

* 

IF  (AT.LE.O.  )  THEN 

MAR05560 
•  MAR05570 

AT=0 

MAR05580 

GOTO  750 

MAR05590 

END  IF 

MAR05600 

IF  (ABS(AT-ATI).GT. .0001)  GOTO  730 

MAR05610 

** 

COMPUTE  THE  EMPIRICAL  BAYES  ATTRITION  RATES 

MAR05620 

750 

DO  770  IM=1 ,NMOS 

MAR05630 

DO  780  IL=1,NL0S 

MAR05640 

DO  790  IG=1 ,NGRD 

MAR05650 

IF  (D(IM,IL,IG).EQ.  1)  THEN 

MAR05660 

DENT  =  AT  +  VT( IM,IL,IG) 

MAR05670 

XTEB( IM, IL, IG)  =  (AT/DENT)*XTB( IM, IL, IG)  + 

MAR05680 

*  ( ( DENT-AT) /DENT)*ZTB 

MAR05690 

ENDIF 

MAR05700 

790 

CONTINUE 

MAR05710 

780 

CONTINUE 

MAR05720 

770 

CONTINUE 

MAR05730 

** 

COMPUTE  THE  EMPIRICAL  BAYES  SHRINKAGE  RATES 

MAR05740 

DO  800  IM=1,NM0S 

MAR05750 

DO  801  IL=1 ,NL0S 

MAR05760 

DO  802  IG=1,NGRD 

MAR05770 

SHRTEBC IM, IL, IG)  =  0 

MAR05780 

IF  (D( IM, IL,IG).  EQ.  1)  THEN 

MAR05790 

SHRTEBC IM.IL.IG)  =  VT( IM,IL,IG)/(AT  +  VT(IM,IL,IG)) 

MAR05800 

ENDIF 

MAR05810 

802 

CONTINUE 

MAR05820 

801 

CONTINUE 

MAR05830 

800 

•VJUUL 

CONTINUE 

MAR03840 

MAR05850 

TC7V7T  ft 
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** 

friWnV 

JAMES -STEIN  SHRINKAGE  RATES.  ** 

MAR05860 

MAR05870 

SST=0 

MAR05880 

SSE=0 

MAR05890 

SSB=0 

MAR05900 

iHr 

COMPUTE  THE  SUMS  OF  SQUARES  (SST,  SSE  AND  SSB) 

MAR059 10 

DO  90  IM=1 ,NMOS 

MAR05920 

DO  95  IL=l,NLOS 

MAR05930 

DO  110  IG=1,NGRD 

MAR05940 

TEMP 1=0 

MAR05950 

TEMP2=0 

MAR05960 

DO  120  IT=1 ,T 

MAR05970 

TEMP 1=TEMP l+( X( IM, IL, IG, IT) -XBB)**2 

MAR05980 

TEMP2=TEMP2+( X( IM , IL, IG, IT) -XB( IM, IL, IG) )**2 

MAR05990 

120 

CONTINUE 

MAR06000 

SST=SST+D( IM , IL , IG) *TEMP 1 

MAR06010 

SSE=SSE+D( IM, IL, IG)*TEMP2 

MAR06020 

110 

CONTINUE 

MAR06030 

95 

CONTINUE 

MAR06040 

90 

CONTINUE 

MAR06050 

SSB=SST-SSE  " 

MAR06060 

★★ 

COMPUTE  THE  JAMES-STEIN  SHRINKAGE  RATES 

MAR06070 

SHRX=(SSE/SSB)*  (KK-3)/(KK*(T-l)+2) 

.  * 

MAR06080 

61 


SHR-DMIN 1 ( SHRX ,  l.DO) 

SHR=DMAX1(SHR  ,  O.DO) 

SKR1=T*SHR 

*************  ******* A******** *********  ***************  ****************** 

**  LIMITED  TRANSLATION  OPTION  &  JAMES -STEIN  ** 

**  ON  TRANSFORMED  SCALE.  ** 

*********  *  ****  *******  ****  *****************  **"*  A  A  A  A  *******  *******  A  ****** 

NN=KK*(T-1) 

DO  125  IM»l,NMOS 
DO  130  IL=1 ,NLOS 
DO  140  IG=1,NGRD 
XJ(IM,IL,IG)=o 
XJ1(IM,IL,IG)=0 
IF  (D(IM,IL,IG).EQ. 1)  THEN 
XJ( IM, IL, IG)=XBB+( 1-SHR)*(XB( IM, IL,IG)-XBB) 
XJ1(IM,IL,IG)=XBB+(1-SHR1)*(XB(IM,IL,IG)-XBB) 

C  ***  XJLT(IM,IL,IG)=XBB+(1-SHR*RH0(XB(IM,IL,IG),XBB,SSE, 

C  ***  SSB,DEE,KK,T,NN))*(XB(IM,IL,IG)-XBB) 

XMLE( IM, IL, IG)=SQRT( .  5+ACINVC IM, IL, IG) ) 

*  *ASIN( 2*PMLE( IM, IL, IG) -1) 

ENDIF 

140  CONTINUE 
130  CONTINUE 
125  CONTINUE 

******************  ***  ***********  **********  *******  ****  **•*•*  ************* 

**  VALIDATION.  RISK  EVALUATED  IN  TRANSFORM  SPACE.  ** 

********************************************************************** 

DO  150  IT=1 ,  NVYR 
RTSCAC IT) =0 
RJ(IT)=0 
RJ1(IT)=0 
RXEB( IT)=0 
RXTEB(  IT)=0 
C  ***  RJLT( IT)=0 
RML( IT)=0 
DO  160  IM= 1 , NMOS 
DO  170  IL=1,NL0S 
DO  180  IG=1,NGRD 

IF(D(IM,IL,IG).  EQ.  1  .AND.  CINV(IM,IL,IG,IT+T).  NE.  0)  THEN 
RTSCAC IT)=RTSCA( IT)+( X( IM , IL , IG ,T+IT) -  XB( IM , IL, IG) )**2 
RJ(IT)=  RJ(IT)+(X(IM,IL,IG,T+IT)-  XJ( IM, IL, IG) )**2 
RJ1(IT)=  RJ1(IT)+(X( IM,IL,IG,T+IT)-  XJ1C IM, IL, IG) )**2 
RXEB( IT)=RXE8( IT)+( X( IM, IL, IG.T+IT) -  XEB( IM, IL, IG) )**2 
RXTEBC IT)=RXTEB( IT)+((X( IM, IL, IG,T+IT)/SQRT(0. 5  + 

CINV( IM , IL , IG , IT+T) ) -  XTEBC IM , IL, IG) )**2) 

/VT( IM, IL, IG) 

RJLT  ( IT)=  RJLT( IT)+(X( IM,IL, IG,T+IT)-XJLT( IM, IL, IG) )**2 
RML  (IT)=  RML  (IT)+(X(IM,IL,IG,T+IT)-XMLE(IM,IL,IG))**2 
ENDIF 
CONTINUE 
CONTINUE 
160  CONTINUE 

RTSCAC IT)=RTSCA( IT)/KV( IT) 

RJ  (IT)*RJ  ( IT)/KV( IT) 

RJ1  (IT)=RJ1  ( IT)/KV( IT) 

RXEBC IT)=RXEB( IT)/KV( IT) 


* 

* 

C  *** 


180 

170 
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MAR06140 
MAR06150 
MAR06160 
MAR06170 
MAR06180 
MAR06190 
MAR06200 
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MAR06610 
MAR06620 
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RXTEB( IT)=RXTEB( IT) /( KV( IT)*T) 

C  ***  RJLT( IT)=RJLT( IT) /KV( IT) 

RML  ( IT)=RML  (IT)/KV(IT) 

RTJ( IT)=RTSCA( IT) /RJ( IT) 

C  ***  RTJLTC IT)=RTSCA( IT) /RJLT( IT) 

C  ***  RSL( IT)=(RJ( IT) -RTSCAC IT) )/ (RJLT( IT) -RTSCA( IT) ) 

150  CONTINUE 

***************  ********  ****************  **  ***********  A**  *************** 

**  DO  RISK  ANALYSIS  ON  THE  ORIGINAL  SCALE.  ** 

******  *********  ***  *****  A*  *************** *****  ********  *************** 

DO  200  IM»1,NM0S 
DO  210  IL=1 ,NLOS 
DO  220  IG=1,NGRD 

**  AVERAGE  CENTRAL  INVENTORY  OVER  ESTIMATION  SET 

CINVS=0 
DO  230  IT-1, T 

AVCINV=AVCINV  +  CINV(IM,IL,IG,IT) 

230  CONTINUE 

AVCINV=AVCINV/T 

TSCA(IM,IL,IG)=D(IM,IL,IG)*SCINV(  XB( IM, IL.IG) ,  AVCINV) 
PJ(IM,IL,IG)=D(IM,IL,IG)*SCINV(  XJ( IM, IL, IG) ,  AVCINV) 
PJKIM.IL, IG)=D(IM,IL,IG)*SCINV(  XJ1( IM,IL, IG)  ,  AVCINV) 
PEB(IM,IL,IG)=D(IM,IL,IG)*SCINV(  XEB( IM.IL.IG) ,  AVCINV) 
PTEB(IM,IL,IG)=D( IM,IL,IG)*Q. 5*( 1.  +  SIN(XTEB(IM,IL,IG))) 

C  ***  PJLT( IM, IL, IG)®D( IM, IL, IG)*SCINV(XJLT( IM.IL.IG) ,  AVCINV) 

IF  (D(IM,IL,IG).EQ. 1)  THEN 
TSCAM( IM.IL.IG) =MILLER(  XB(IM,IL,IG) ,  AVCINV  ) 

PJM( IM.IL.IG )=MILLER(  XJ( IM.IL.IG) ,  AVCINV  ) 

PJ1MC IM, IL, IG)=MILLER(  XJ1( IM.IL.IG) ,  AVCINV  ) 

PEBM( IM.IL.IG )=MILLER(  XEB(IM, IL.IG) ,  AVCINV  ) 

XMIL  =  XTEBC IM.IL.IG)  *  SQRTfO. 5  +  AVCINV) 

PTEBMC IM, IL, IG)«MILL£R(  XMIL  ,  AVCINV  ) 

END  IF 

220  CONTINUE 
210  CONTINUE 
200  CONTINUE 

*********  * ** ********  ******  **  *  **  ***********  ******  ** **  A  *  A  *****  ********** 

**  VALIDATION.  RISK  EVALUATED  IN  ORIGINAL  SPACE.  ** 

**  ************************ ** **** **** *  *************** * ******** *★* *****  * 
DO  240  IT=1,  NVYR 
RRTSCA( IT)=0 
RRTSAMC IT)=0 
RRJ( IT)=0 
RRJM( IT)=0 
RRPEB( IT)=0 
RRPTEB(IT)=0 
RRPEBM(IT)=0 
RRPTBMC IT)=0 
RRJ1(IT)=0 
RRJ1M( IT)=0 
C  ***  RRJLTC IT)=0 
RRML( IT)=0 
KP=0 

DO  250  IM=1,NM0S 
DO  260  IL-l.NLOS 
DO  270  IG-l.NGRD 


MAR06650 
MAR06660 
MAR06670 
MAR06680 
‘  MAR06690 
MAR06700 
MAR06710 
MAR06720 
MAR06730 
MAR06740 
MAR06750 
MAR06760 
MAR06770 
MAR06780 
MAR06790 
MAR06800 
MAR06810 
MAR06820 
MAR06830 
MAR06840 
MAR06850 
MAR06860 
MAR06870 
MAR06880 
MAR06890 
MAR06900 
MAR06910 
MAR06920 
MAR06930 
MAR06940 
MAR06950 
MAR06960- 
MAR06970 
MAR06980 
MAR06990 
MAR07000 
MAR07010 
MAR07020 
MAR07030 
MAR07040 
MAR07050 
MAR07060 
MAR07070 
MAR07080 
MAR07090 
MAR07100 
MAR07110 
MAR07120 
MAR07130 
MAR07140 
MAR07150 
MAR07160 
MAR07170 
MAR07180 
MAR07190 
MAR07200 


-  * 

63 


* 


* 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

c  *** 
c  * 
** 

* 

* 


CC*C INV( IM , IL , IG , T+IT) 

IF(CC  .NE. 0)  THEN 

AR=Y( IM.IL.IG, T+IT)/CC 


ELSE 


AR=0 


END  IF 

IF(D( IM, IL, IG)  .NE.  0  .AND.  PMLE(IM,IL,IG).  NE. 0 
.AND.  PMLE( IM,IL,IG).  NE.  1)  THEN 


KP=KP+1 

RRML( IT)=RRML( IT)+( CC*( AR-PMLE( IM, IL, IG) )**2) 

/(PMLE( IM,IL,IG)*( 1.  -  PMLE(IM,IL,IG))) 


ENDIF 

PRODUCE  CHI-SQUARE  MEASURE  OF  EFFECTIVENESS 

RRTSCAC IT)*  RRTSCA( IT)+CHISQR(D( IM, IL, IG) ,TSCA( IM, IL, IG) , 

CC,AR) 

RRTSAM( IT)*  RRTSAMC IT)+CHISQR(D( IM, IL, IG) ,TSCAM( IM, IL, IG) , 

CC,AR) 

RRJ  (IT)*  RRJ  ( IT)+CHISQR(D( IM,IL,IG) ,PJ  (IM,IL,IG), 

CC.AR) 

RRJM  (IT)*  RRJM  (IT)+CHISQR(D( IM.IL.IG) ,PJM  (IM,IL,IG) , 

CC.AR) 

RRJ1  (IT)*  RRJ1  (IT)+CHISQR(D(IM,IL,IG),PJ1  (IM.IL.IG), 

CC.AR) 

RRJ1M  ( IT)3  RRJ1M  ( IT)+CHISQR(D( IM, IL,IG) ,PJ1M( IM.IL.IG) , 

CC,AR) 

RRPEB  (IT)*  RRPEB  (IT)+CHISQR(D( IM.IL.IG) ,PEB  (IM.IL.IG), 

CC.AR) 

RRPTEB  (IT)*  RRPTEB  (IT)+CHISQR(D( IM.IL.IG) ,PTEB  (IM.IL.IG), 

CC.AR) 

RRPEBM( IT)*  RRPEBM( IT)+CHISQR(D( IM, IL, IG) ,PEBM( IM, IL, IG) , 

CC.AR) 

RRPTBM( IT)=  RRPTBM( IT)+CHISQR(D( IM , IL, IG) ,PTEBM( IM, IL, IG) , 

CC.AR) 

RRJLT  (IT)*  RRJLT  ( IT)+CHISQR(D( IM, IL, IG) ,PJLT(IM, IL, IG) , 

CC.AR) 

PRODUCE  MEAN  ABSOLUTE  DEVIATION  MEASURE  OF  EFFECTIVENESS 

IF(D( IM.IL.IG)  .NE.  0  .AND.  PMLE( IM, IL, IG).  NE.  0 
.AND.  PMLE(IM,IL,IG). NE. 1)  THEN 
CALL  UOM(  CC,  AR,  PMLE( IM.IL.IG) ,  MOPMLE(IT), 

MUPMLE(IT),  KVOPML(IT),  KVUPML(IT),  KVDPML(IT)) 


ENDIF 


IF(CC. NE. 0. )  THEN 

CALL  UOM(  CC,  AR,  TSCA( IM.IL.IG) ,  MOTSCA(IT), 

*  MUTSCA(IT),  KVOTSA(IT),  KVUTSA(IT),  KVDTSA(IT)) 

CALL  UOM(  CC,  AR,  TSCAM( IM.IL.IG) ,  MOTSAM(IT), 

*  MUTSAM(IT),  KVOTSM(IT),  KVUTSM(IT),  KVDTSM(IT)) 

CALL  UOM(  CC,  AR,  PJ( IM, IL, IG) ,  MOPJ(IT), 

*  MUPJ(IT),  KVOPJ(IT),  KVUPJ(IT),  KVDPJ(IT)) 

CALL  UOM(  CC,  AR,  PJM( IM, IL, IG) ,  MOPJM(IT), 

*  •  MUPJM(IT),  KVOPJM(IT),  KVUPJM(IT),  KVDPJM(IT)) 

CALL  UOM(  CC,  AR,  PJ1( IM, IL, IG) ,  MOPJl(IT), 

G4 
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MAR07580 
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MUPJl(IT),  KVOPJl(IT),  KVUPJl(IT),  KVDPJl(IT)) 

CALL  UOM(  CC,  AR,  PJ1M( IM, IL, IG) ,  MOPJIM(IT), 

MUPJIM(IT),  KVOPIM(IT),  KVUPIM(IT),  KVDPIM(IT)) 


CALL  UOM(  CC,  AR,  PEB( IM, IL, IG) ,  MOPEB(IT), 

MUPEB(IT),  KVOPEB(IT),  KVUPEB(IT),  KVDPEB(IT)) 

CALL  UOM(  CC,  AR,  PTEB( IM, IL, IG) ,  MOPTEB(IT), 

MUPTEB(IT),  KVOPTB(IT),  KVUPTB(IT),  KVDPTB(IT)) 


CALL  UOM(  CC,  AR,  PEBM( IM, IL, IG) ,  MOPEBM(IT), 

*  MUPEBM(IT),  KVOPBM(IT),  KVUPBM(IT),  KVDPBM(IT)) 

CALL  UOM(  CC,  AR,  PTEBM( IM,IL,IG) ,  MOPTBM(IT), 

*  MUPTBM(IT),  KVOBTM(IT),  KVUBTM( IT) ,  KVDBTM(IT)) 
ENDIF 

270  CONTINUE 
260  CONTINUE 

250  CONTINUE 

IF  (KP  .EQ.  0)  THEN 
FF  =  0 

ELSE 

FF=FLOAT( KV( IT) ) / FLOAT ( KP ) 

ENDIF 

RRML( IT)=RRML( IT)*FF 
240  CONTINUE 

**  COMPUTE  TIME  AVERAGE  MAD 

**  MEASURE  OF  EFFECTIVENESS 

DO  300  IT=1,  NVYR 

MUPMLE(IT)  *  MUPMLE( IT)/KV( IT) 

MUTSCA(IT)  *  MUTSCA( IT) /KV( IT) 

MUTSAM(IT)  =  MUTSAM( IT) /KV( IT) 

MUPJ  (IT)  =  MUPJ  ( IT)/KV( IT) 

MUPJM  (IT)  =  MUPJM( IT) /KV( IT) 

MUPJ1  (IT)  -  MUPJ1( IT) /KV( IT) 

MUPJIM(IT)  =  MUPJ1M( IT)/KV( IT) 

MUPEB  (IT)  =  MUPEB( IT)/KV( IT) 

MUPTEB  (IT)  =  MUPTEB( IT)/KV( IT) 

MUPEBM(IT)  *  MUPEBM( IT)/KV( IT) 

MUPTBM(IT)  *  MUPTBM( IT)/KV( IT) 

MOPMLE(IT)  =  MOPMLE( IT) /KV( IT) 

MOTSCA(IT)  »  MOTSCA( IT) /KV( IT) 

MOTSAM(IT)  =  MOTSAM( IT) /KV( IT) 

MOPJ  (IT)  *  MOPJ  ( IT)/KV( IT) 

MOPJM  (IT)  =  MOPJM( IT) /KV( IT) 

M0PJ1  (IT)  =  M0PJ1( IT) /KV( IT) 

MOPJIM(IT)  *  M0PJ1M( IT) /KV( IT) 

MOPEB  (IT)  =  MOPEB( IT) /KV( IT) 

MOPTEB  (IT)  =  MOPTEB( IT)/KV( IT) 

MOPEBM(IT)  *  MOPEBM( IT)/KV( IT) 

MOPTBM(IT)  ■  MOPTBM( IT) /KV( IT) 

MDPMLE(IT)  -  MUPMLE(IT)  +  MOPMLE(IT) 

MDTSCA(IT)  =  MUTSCA(IT)  +  MOTSCA(IT) 

MDTSAM(IT)  *  MUTSAM(IT)  +  MOTSAM(IT) 

65 


MAR07770 

MAR07780 

MAR07790 

MAR07800 

MAR07810 

MAR07820 

MAR07830 

MAR07840 

MAR07850 

MAR07860 

MAR07870 

MAR07880 

MAR07890 

MAR07900 

MAR07910 

MAR07920 

MAR07930 

MAR07940 

MAR07950 

MAR07960 

MAR07970 

MAR07980 

MAR07990 

MAR08000 

MAR08010 

MAR08020 

MAR08030 

MAR08040 

MAR08050 

MAR08060 

MAR08070 

MAR08080 

MAR08090 

MAR08100 

MAR08110 

MAR08120 

MAR08130 

MAR08140 

MAR08150 

MAR08160 

MAR08170 

MAR08180 

MAR08190 

MAR08200 

MAR08210 

MAR08220 

MAR08230 

MAR08240 

MAR08250 

MAR08260 

MAR08270 

MAR08280 

MAR08290 

MAR08300 

MAR08310 

MAR08320 


MDPJ  (IT)  *  MUPJ  (IT)  +  MOPJ  (IT) 

MDPJM  (IT)  *  MUPJM  (IT)  +  MOPJM  (IT) 

MDPJ1  (IT)  *  MUPJ1  (IT)  +  M0PJ1  (IT) 

MDPJ1M( IT)  »  MUPJ1M( IT)  +  M0PJ1M( IT) 

MDPEB  (IT)  »  MUPEB  (IT)  +  MOPEB  (IT) 

MDPTEB  (IT)  =  MUPTEB  (IT)  +  MOPTEB  (IT) 
MDPEBM(IT)  *  MUPEBM(IT)  +  MOPEBM(IT) 
MPPTBM(IT)  *  MUPTBM(IT)  +  MOPTBM(IT) 

300  CONTINUE 

*********************************************** 

**  PRINT  SECTION  FOR  ALL  OUTPUT. 

*******  *  *  **** **  *  *  *  *  ********  *  *******  *  *  *  *****  ** ** 


**  PRINT  GROUP  CHARACTERISTICS 

VRITE(6, ' ( IX, ' ' TEN  OR  LESS  BAYES  VARIANCE  ITERATIONS  "  ) ' ) 
WRITER,’ (IX  ,12,"  BAYES  VARIANCE  ITERATION  S)  ")* )  IC 
VRITE(  6 , '  (  IX,  ’ '  XBB  ",F12.  6)*)  XBB 
WRITE(6.905)  T,KK,SSB,SSE,SHRX.SHR,ZB,ZTB,PZTB,A,AT 

905  FORMATS  GROUP  CHARACTERISTICS' ,  /,'  T=' ,16,  / ,  K=' ,16, 

*  /,'  SSB=' ,F12. 4,  /,'  SSE=',F12.4,  /,’  SHRXa’ ,F12. 4, 

*  /,'  SHR=',F12.4,  /,'  ZB=',F12.4,  /,’  ZTBa',F12.4, 

*  /,’  PZTB=’ ,F12. 4,  /,'  A=',F12. 4,  /,'  ATa',F12.4) 

**  PRINT  RESULTS  FOR  TRANSFORM  SPACE 

WRITE(6,*)  '  RESULTS  FOR  TRANSFORM  SPACE' 

WRITE(6,906)  (IVYR(I),  1=1,  NVYR) 

906  F0RMAT(14X, 10112) 

WRITE(6.907)  ' RTSCA  * ,  (RTSCA(IT),  IT=1,NVYR) 

907  FORMAT('  ' ,A14, 10F12. 4) 

WRITE(6,907)  'RJ  *,  (RJ  (IT),  IT=1,NVYR) 

WRITE(6,907)  'RJl  ',  (RJ1  (IT),  IT=1,NVYR) 

VRITE(6 ,907)  'RXEB  '  (RXEB  (IT),  IT=1,NVYR) 

WRITE(6,907)  ' RXTEB  f,  (RXTEB  (IT),  IT=1,NVYR) 

C  ***  WRITE(6,907)  'RJLT  ’,  (RJLT  (IT),  IT=1,NVYR) 

WRITE(6,907)  'RML  ',  (RML  (IT),  IT»1,NVYR) 

C  ***  WRITE(6,907)  'RSL  ’,  (RSL  (IT),  IT»1,NVYR) 

C  ***  WRITE(6,907)  'RTJ  ',  (RTJ  (IT),  IT=1,NVYR) 

C  ***  WRITE(6,907)  ' RTJLT  ',  (RTJLT(IT),  IT=1,NVYR) 

**  PRIOT  RESULTS  FOR  ORIGINAL  SPACE 


WRITE(6,*)  '  ' 

WRITE(6,*)  ’  RESULTS  FOR  ORIGINAL  SPACE' 


WRITE(6,907)  ' RRML 


(RRML  (IT),  IT=1 ,NVYR) 


WRITE(6,907)  ' RRTSCA'  , 
WRITE( 6,907)  'RRTSAM' , 
WRITE(6,907)  'RRJ  ', 
WRITE(6 ,907)  'RRJM  ', 
WRITE(6,907)  'RRJl  ', 
WRITE(6,907)  'RRJ1M  ', 
WRITE(6,907)  'RRPEB 
WRITE( 6,907)  ’ RRPTEB  f, 
WRITE(6,907)  'RRPEBM' , 
WRITE(6,907)  'RRPTBM' , 
C  ***  WRITE(6,907)  'RRJLT  ', 
WRITE(6.908)  *KV  \ 
908  FORMAT(  ',A14. 10112) 
¥RITE(6, ' (9X,  KP  " 
WRITE(6,*)  '  ' 


(RRTSCA(IT),  IT=1,NVYR) 
(RRTSAM(IT),  IT=1,NVYR) 
(RRJ  (IT),  IT=1 ,NVYR) 
(RRJM  (IT),  IT=1,NVYR) 
(RRJl  (IT),  IT=1 ,NVYR) 
(RRJ1M  (IT),  IT=1 ,NVYR) 
(RRPEB(IT),  IT=1,NVYR) 
(RRPTEB(IT),  IT=1,NVYR) 
(RRPEBM(IT),  IT=1 ,NVYR) 
(RRPTBM(IT),  IT=1,NVYR) 
(RRJLT(IT),  IT=1,NVYR) 
(KV  (IT),  IT=1 ,NVYR) 

,112)')  KP 


PRINT  MAD  CALCULATIONS 


MAR08330 
MAR06340 
MAR08350 
MAR08360 
•  MAR08370 
MAR08380 
MAR08390 
MAR08400 
MAR08410 
MAR08420 
MAR08430 
MAR08440 
MAR08450 
MAR08460 
MAR08470 
MAR08480 
MAR08490 
MAR08500 
MAR08510 
MAR08520 
MAR08530 
MAR08540 
MAR08550 
MAR08560 
MAR08570 
MAR08580 
MAR08590 
MAR08600 
MAR08610 
MAR08620 
MAR08630 
MAR08640 
MAR08650 
MAR08660 
MAR08670 
MAR08680 
MAR08690 
MAR08700 
MAR08710 
MAR08720 
MAR08730 
MAR08740 
MAR08750 
MAR08760 
MAR08770 
MAR08780 
MAR08790 
MAR08800 
MAR08810 
MAR08820 
MAR08830 
MAR08840 
MAR08850 
MAR08860 
MAR08870 
MAR08880 


WRITE ( 6 , 1 ( 2X  1 1 PMLE 1 1 ) 1 ) 

CALL  WRMAD ( I  VYR , NVYR . MUPMLE , KVUPML , MOPMLE , KVOPML , MDPMLE , KVDPML) 
WRITEC 6 , ' ( 2X,  " TSCA' ' ) ' ) 

CALL  WRMAD ( I VYR , NVYR . MUTSCA , KVUTSA , MOTSCA , KVOTSA  ,MDTSCA , KVDTSA) 
WRITE(  6 , '  (  2X ,  "  TSCAM 1 ' ) '  ) 

CALL  WRMAD ( I VYR , NVYR . MUTS AM , KVUTSM , MOTSAM , KVOTSM , MDTSAM , KVDTSM) 
WRITE(6 , ' ( 2X, ' 1 PJ' ' ) '  ) 

CALL  WRMAD ( I VYR , NVYR , MUP J , KVUPJ , MOP J , KVOP J , MDP J , KVDP J ) 
WRITE(6,'(2X,"PJM")') 

CALL  WRMAD ( I VYR , NVYR , MUPJM , KVUP JM , MOPJM, KVOP JM , MDP JM , KVDP JM) 
WRITEC  6 » ' ( 2X, ' ' PJ1 "  ) ' ) 

CALL  WRMADC I VYR, NVYR, MUP J1 .KVUPJl .MOPJl .KVOPJl ,MDPJ1 .KVDPJl) 
WRITEC  6,' (2X,"PJlM,f)') 

CALL  WRMADC I VYR , NVYR , MUP JIM , KVUP 1M , MOP J 1M , KVOP 1M , MDP J 1M , KVDP1M) 
WRITEC  6 , ' C  2X , ' ' PEB 1 1 ) 1 ) 

CALL  WRMADC I VYR , NVYR . MUPEB , KVUPEB , MOPEB , KVOPE8 , MDPEB , KVDPEB) 
WRITEC 6 C 2X,  ”  PEBM' f )' ) 

CALL  WRMADC I VYR , NVYR , MUPEBM , KVUPBM , MOPEBM , KVOPBM , MDPEBM , KVDPBM) 
WRITEC 6, ' ( 2X, ' ' PTEB' ' ) '  ) 

CALL  WRMADC I VYR , NVYR . MUPTEB , KVUPTB , MOPTEB , KVOPTB .MDPTEB , KVDPTB ) 
WRITEC 6,' C2X,  "  PTEBMf' )' ) 

CALL  WRMADC I VYR . NVYR , MUPTBM , KVUBTM , MOPTBM , KVOBTM , MDPTBM , XVDBTM) 
WRITEC6,*)  '  ' 

**  PRINT  THE  CENTRAL  INVENTORY  MATRIX 

WRITEC 6,*)  ' CINV  MATRIX' 

DO  1006  IT=1 ,NYR 
WRITEC 6,*)  YEAR', IT 
DO  1005  IL=1 ,NLOS 

WRITEC 6 , 1007 ) C CINVC IM , IL, 1 , IT) , IM=1 , NMOS ) 

1005  CONTINUE 

1006  CONTINUE 

1007  F0RMATC15C1X.F7. 2)) 

WRITEC 6,*)  '  ' 

**  PRINT  THE  ATTRITION  MATRIX 

WRITEC 6,*)  'Y  MATRIX' 

DO  1016  IT=1 ,NYR 
WRITEC6,*)  ‘YEAR’, IT 
DO  1015  IL=1,NL0S 

WRITEC  6 , 1007 ) C  YC IM , IL , 1 , IT) , IM=1 , NMOS ) 

1015  CONTINUE 

1016  CONTINUE 

WRITEC 6,*)  '  ' 

**  PRINT  ALL  ARRAYS 

WRITEC6,'C2X,  "ARRAY  TSCA")') 

CALL  ARRAYC  MXMOS , MXLOS , MXGRD , NMOS , NLOS , NGRD , TSCA ) 

WRITEC6,'C2X,  "ARRAY  TSCAM")’) 

CALL  ARRAYC MXMOS, MXLOS, MXGRD, NMOS, NLOS, NGRD, TSCAM) 

WRITEC6,'C2X,  "ARRAY  SHREB")') 

CALL  ARRAYC  MXMOS , MXLOS , MXGRD , NMOS , NLOS , NGRD , SHREB ) 

WRITEC 6,' C2X," ARRAY  SHRTEB "  ) ' ) 

CALL  ARRAYC  MXMOS , MXLOS , MXGRD , NMOS , NLOS , NGRD , SHRTEB ) 
WRITEC6,'C2X,  "ARRAY  PMLE "  ) '  ) 

CALL  ARRAYC MXMOS, MXLOS, MXGRD, NMOS, NLOS, NGRD, PMLE) 

WRITEC  6, 'C2X,M  ARRAY  PJ")') 

CALL  ARRAYC  MXMOS , MXLOS , MXGRD , NMOS , NLOS , NGRD , P J) 

WRITEC6,'C2X,  "ARRAY  PJM")') 

-'■T 
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MAR08890 

MAR08900 

MAR08910 

MAR08920 

MAR08930 

MAR08940 

MAR08950 

MAR08960 

MAR08970 

MAR08980 

MAR08990 

MAR09000 

MAR09010 

MAR09020 

MAR09030 

MAR09040 

MAR09050 

MAR09060 

MAR09070 

MAR09080 

MAR09090 

MAR09100 

MAR09110 

MAR09120 

MAR09130 

MAR09140 

MAR09150 

MAR09160 

MAR09170 

MAR09180 

MAR09190 

MAR09200 

MAR09210 

MAR09220 

MAR09230 

MAR09240 

MAR09250 

MAR09260 

MAR09270 

MAR09280 

MAR09290 

MAR09300 

MAR09310 

MAR09320 

MAR09330 

MAR09340 

MAR09350 

MAR09360 

MAR09370 

MAR09380 

MAR09390 

MAR09400 

MAR09410 

MAR09420 

MAR09430 

MAR09440 


CALL  ARRAY ( MXMOS , MXLOS , MXGRD , NMOS , NLOS , NGRD , PJM ) 
WRITE(6,'(2X,  "ARRAY  PJl")') 

CALL  ARRAY ( MXMOS , MXLOS , MXGRD . NMOS , NLOS , NGRD , P J1 ) 
WRITE(  6 , '  (  2X ,  "  ARRAY  PJlM")f) 

CALL  ARRAY ( MXMOS , MXLOS , MXGRD , NMOS , NLOS , NGRD , P J 1M) 
WRITE( 6 , ' ( 2X ,  " ARRAY  PEB1 ' )') 

CALL  ARRAY ( MXMOS , MXLOS , MXGRD . NMOS , NLOS , NGRD , PEB ) 
WRITE(6,'(2X,  "ARRAY  PTEB")f) 

CALL  ARRAY ( MXMOS , MXLOS , MXGRD .NMOS, NLOS , NGRD , PTEB ) 
WRITE(6, ' (2X, ' 'ARRAY  PEBM' '  ) '  ) 

CALL  ARRAY( MXMOS , MXLOS , MXGRD , NMOS , NLOS , NGRD , PEBM) 
WRITE(6,'(2X, "ARRAY  PTEBM")’) 

CALL  ARRAY ( MXMOS , MXLOS . MXGRD , NMOS , NLOS , NGRD , PTEBM) 
WRITE( 6 , ' ( 2X , 1 ' ARRAY  V1')’) 

CALL  ARRAY ( MXMOS , MXLOS , MXGRD , NMOS , NLOS , NGRD , V) 
WRITE(6,'(2X,  "ARRAY  XB")') 

CALL  ARRAY( MXMOS, MXLOS, MXGRD, NMOS, NLOS, NGRD, XB) 
WRITE(  6 , '  (  2X,  "  ARRAY  VT")') 

CALL  ARRAY( MXMOS , MXLOS , MXGRD , NMOS , NLOS , NGRD , VT) 
WRITE(  6 , '  (  2X , '  1  ARRAY  XTB")') 

CALL  ARRAY ( MXMOS , MXLOS , MXGRD , NMOS , NLOS , NGRD , XTB ) 
WRITE(  6 , '  (  2X ,  "  ARRAY  XJ"  ) '  ) 

CALL  ARRAY(MXMOS, MXLOS, MXGRD, NMOS, NLOS, NGRD, XJ) 
WRITE(  6 , '  (  2X,  "  ARRAY  XJl")') 

CALL  ARRAY( MXMOS , MXLOS , MXGRD , NMOS , NLOS , NGRD , XJl ) 
WRITE(  6 , '  (  2X ,  "  ARRAY  XEB")') 

CALL  ARRAY ( MXMOS , MXLOS , MXGRD , NMOS , NLOS , NGRD , XEB ) 
WRITE(  6 , '  (  2X,  "  ARRAY  XTEB")') 

CALL  ARRAY ( MXMOS , MXLOS , MXGRD , NMOS , NLOS , NGRD , XTEB ) 


SUBROUTINE  WRMAD( IYR,NYR,MU,KU,MO,KO,MD,KD) 

***************  it*******#**-****** Hrirtrh irtrlrfc  /rfr******  ft  ******  A  ft  ft  *  ft**** 

INTEGER  IYR(NYR) ,KU(NYR) ,KO(NYR) ,KD(NYR) 

REAL  MU(NYR) ,MO(NYR) ,MD(NYR) 

WRITE(6, ' (/9X, ' 'DEVIATION' ' ,T26 ,I2,T42 ,I2,T59 , 12/) ' ) 

*  (IYR(IW) , IW-l.NYR) 

VRTTEC6, '(9X,  "UNDERAGE1  '  ,T22,F9.  6 '  (  "  ,12 ,'  '  )  "  ,T38,F9.  6, 

*  "  (  "  ,12,  "  )"  ,T55  ,F9.  6,  '  (  "  ,12, f  ’  )  "  )  f ) 

*  MU(1) ,KU( 1) ,MU(2) ,KU(2) ,MU(3) ,KU(3) 

WRITE(6,'(9X,  "OVERAGE1  ’  ,T22,F9.6,"  (  "  ,12,")"  ,T38,F9.6, 

*  "  (",I2,")",T55,F9.6,"  (", 12, ")")') 

*  MO( 1) ,KO( 1) ,MO( 2) ,KQ(2) ,MO( 3) ,KO( 3) 

WRITE(6, '  (9X, '  'TOTAL' '  ,T22,F9.  6,"  ("  .12,")"  .T38,F9.  6, 

*  "  (",I2,")",T55,F9.6,"  (  f '  ,  12,  "  )  f '  )  *  ) 

*  MD( 1) ,KD( 1) ,MD(2) ,KD( 2) ,MD( 3) ,KD( 3) 

WRITE(6, '(/)') 

RETURN 

END 


SUBROUTINE  ARRAY ( MM , ML , MG , NM , NL , NG , OB J 1 ) 


REAL  OBJl(MM,ML,MG) 
DO  10  IG1=1,NG 


MAR09450 
MAR09460 
MAR09470 
MAR09480 
*  MAR09490 
MAR09S00 
MAR09510 
MAR09520 
MAR09530 
MAR09540 
MAR09550 
MAR09560 
MAR09S70 
MAR09580 
MAR09590 
MAR09600 
MAR09610 
MAR09620 
MAR09630 
MAR09640 
MAR09650 
MAR09660 
MAR09670 
MAR09680 
MAR09690 
MAR09700 
MAR09710 
MAR09720 
MAR09730 
MAR09740 
MAR09750 
MAR09760 
MAR09770 
MAR09780 
MAR09790 
MAR09800 
MAR09810 
MAR09820 
MAR09830 
MAR09840 
MAR09850 
MAR09860 
MAR09870 
MAR09880 
MAR09890 
MAR09900 
MAR09910 
MAR09920 
MAR09930 
MAR09940 
'  MAR09950 

MAR09960 
r  MAR09970 
MAR09980 
MAR09990 


DO  10  ILl-l.NL 

10  WRITE(6 . f  (2X,  15(F7. 3 , IX) ) ' )  (0BJ1( IM1,IL1,IG1) ,IM1=1,NM) 

WR1TE(6, '(/)') 

RETURN 

END 

A*******1***********  ****************"**  A**  A*  *****************  *********** 

SUBROUTINE  FTT(CINV,Y,  X) 

********************************************************************** 
REAL  CINV ,X,Y 
TEMP=Y/( 1+CINV) 

TEMP1=( 1.  +Y)/( 1+CINV) 

TEMP=SQRT(. 5+CINV)*.  5*(ASIN( -1.  +2. *TEMP)  + 

*  ASIN( -1. +2. *TEMP1) ) 

X=TEMP 

END 

********* *********************** *********************** *************** 

SUBROUTINE  UOM( VINV , ACT , OBJ , OVER , UNDER , KO , KU , KD) 
**************************  ******************************************** 

REAL  OBJ 

IF  (ACT  .GT.  OBJ)  THEN 
UNDER  =  UNDER+VINV*( ACT-OBJ) 

KU=KU+1 

ELSE 

OVER  =  OVER  +VINV*( OBJ-ACT) 

KO=KO+l 
END  IF 
KD=KD+1 
RETURN 
END 

*  ****  *  A  **  ***************  ******  *******  *** *  **********  *********  ***  ******* 

REAL  FUNCTION  RHO(XB,XBB,SSE,SSB,DEE,KK,T,NN) 

***************** a **********  ****  *  *  ************************************ 

INTEGER  T 

REAL*8  SIGSQ,  V,  U 
SIGSQ=SSE/(NN+2) 

V=SSB/T 

U=KK*(KK-3)*(XB-XBB)**2/(V*(KX-1)) 

RH0=DMIN1( 1.  DO,  DEE/DSQRT(U)) 

END 

********************  ******  *******  *  A  *  *  ******  ****  ****  A  ****************** 

REAL  FUNCTION  MILLER(  X,  AVGINV  ) 

*****  ******** *******  *  *** *  *  *********  *********  ***  *  ***********  *********** 
TM  =  X/SQRT( AVGINV+. 5)+l. 570796 
TEMPI  =  ASIN(SQRT( 1.  /(AVGINV+1))) 

IF  (TM  .  LT.  TEMPI)  THEN 
MILLER  -  0.  0 
GOTO  1 
ENDIF 

IF  (TM  .GT.  3.  14159-TEMPI)  THEN 
MILLER  *  1.  0 
GOTO  1 
ENDIF 

IF  (COS(TM)  .EQ.  0.0)  THEN 
MILLER  =  0.5 
GOTO  1 
ENDIF 


MAR10000 

MAR10010 

MAR10020 

MAR10030 

MAR10040 

MAR10050 

MAR10060 

MAR10070 

MAR10080 

MAR10090 

MAR10100 

MAR10110 

MAR10120 

MAR10130 

MAR10140 

MAR10150 

MAR10160 

MAR10170 

MAR10180 

MAR10190 

MAR10200 

MAR10210 

MAR10220 

MAR10230 

MAR10240 

MAR10250 

MAR10260 

MAR10270 

MAR10280 

MAR10290 

MAR10300 

MAR10310 

MAR10320 

MAR10330 

MAR10340 

MAR10350 

MAR10360 

MAR10370 

MAR10380 

MAR10390 

MAR10400 

MAR10410 

MAR10420 

MAR10430 

MAR10440 

MAR10450 

MAR10460 

MAR10470 

MAR10480 

MAR10490 

MAR10500 

MAR10510 

MAR10520 

MAR10530 

MAR10540 

MAR10550 
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IF  (COS(TM)  .  GT.  0.0)  SIGNUM  =  1. 

IF  (COS(TM)  .LT.  0.0)  SIGNUM  =  -1. 

TEMP2=1-(SIN(TM)+(SIN(TM)-(1. /SIN(TM)))/AVGINV)**2 
IF  (TEMP2.  LT.  0.  0)  TEMP2=0.  0 
MILLER*.  5*(1.  -SIGNUM*TEMP2**. 5) 

1  END 

********************************************************************** 
REAL  FUNCTION  SCINV(X.A) 

A****************'*'*******'*****'******************'***********'*********** 

**  SCALE  INVERSION  (TRANSFORM  — >  ORIGINAL  SCALE) 

R=X/SQRT(A+.5) 

IF(R  .LT.  -1.570796)  THEN 
SCINV=0. 

ELSE  IF(R  .GT.  1.570796)  THEN 
SCINV=1. 

ELSE 

SCINV=.  5*( 1.  +SIN(R) ) 

END  IF 
END 

*************  **A'**V  v*  ***********  *■**  *  *  *■■***  *******  **************  ******  * 

REAL  FUNCTION  VAR(  ZZ  ) 

****************** **************************************************** 
DATA  AA/1. 683S/,  B1/-.  8934/,  B2/.8991/ 

IF  (  ZZ  .GT.  2.2)  THEN 
VAR  =  1. 

ELSE 

IF  (  ZZ  .LT.  1.001)  THEN 
VAR  =  .05 
ELSE 

VAR  =  AA*(ZZ**B1)*(ZZ-1.  )**B2 
ENDIF 
END  IF 
END 

AAAA**A-*A"******A'AA'*AAA  A***' A  A  **■*****'*******  **'**■*■*  A***  A'**  A  *******  A  ****** 

REAL  FUNCTION  CHISQR(ID,OBJ,INV,ACT) 

*  *  ***  *  ***  ***  *  A  *  A  A  A  A  A  A  A  A  *  **  A  A **  A  >■■******  A**  ************  *******  *  A  A  A  ***  *** 

REAL  OBJ, INV 

IF  (ID.  NE.  0  .AND.  OBJ.  NE.  0.  .AND.  OBJ.  NE.  1.  ) 

*  CHISQR  =  (INV*(ACT-OBJ)**2)/(OBJ*(l. -OBJ)) 

END 

******  A'A  A ,  A ***************************  A  *** A************** A*** ********** 

SUBROUTINE  READER( MOS , LOS , GRADE , YR , INV , Y ,  NMOS.SMOS, 

*  NGRD.SGRD,  SLOS1.SLOS2,  ST1.ST2,  IM,IL,IG,IT,  IEOF) 

**********  **************  *************  ***  AAA  ********  *  **'*'*  ******  A  A  ****** 

**  READ  RECORD  AND  ACCUMULATE  LOSSES 

INTEGER  MOS,  LOS,  GRADE,  YR,  INV,  L0SS(8) 

REAL  Y 

INTEGER  ST1 ,  ST2 
INTEGER  SLOSl,  SLOS2 
INTEGER  SMOS(40) ,  NMOS 
INTEGER  SGRD(IO),  NGRD 
ICNT=0 

1  READ( 1 , 100 ,END=999)  YR, MOS, GRADE, LOS, INV, (LOSS(I) ,1=1,8) 

100  F0RMAT(4I2,9I5) 

**  CHECK  IF  RECORD  MEETS  SELECTION  CRITERIA.  OTHERWISE  REJECT. 

IM=0 


MAR10560 
MAR  105  70 
MAR10560 
MAR10590 
MAR10600 
MAR10610 
MAR10620 
MAR10630 
MAR10640 
MAR10650 
MAR10660 
MAR10670 
MAR10680 
MAR10690 
MAR10700 
MAR10710 
MAR10720 
MAR10730 
MAR10740 
MAR10750 
MAR10760 
MAR10770 
MARIO  780 
MAR10790 
MAR10800 
MAR10810 
MAR10820 
MAR10830 
MAR10840 
MAR10850 
MAR10860 
MAR10870 
MAR10880 
MAR10890 
MAR10900 
MAR1Q910 
MAR10920 
MAR10930 
MAR10940 
MAR10950 
MAR10960 
MAR10970 
MAR10980 
MAR10990 
MAR11000 
MAR11010 
MAR11020 
MAR11030 
MARI 1040 
MAR11050 
MARI 1060 
MAR11070 
MAR11080 
MAR11090 
MARI 1100 
MAR11110 
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DO  10  I=1,NM0S 

MAR11120 

IF(MOS  .EQ.  SMOS(I))  THEN 

MAR11130 

IM=I 

MAR11140 

GO  TO  20 

MARI 1150 

END  IF 

‘MAR11160 

10 

CONTINUE 

MAR11170 

GO  TO  1 

MARI 1180 

★ 

MAR11190 

20 

CONTINUE 

MAR11200 

DO  30  1=1 ,NGRD 

MAR11210 

IF( GRADE  .EQ.  SGRD(I))  THEN 

MAR11220 

IG=I 

MAR11230 

GO  TO  40 

MAR11240 

END  IF 

MAR11250 

30 

CONTINUE 

MAR11260 

GO  TO  1 

MAR11270 

40 

CONTINUE 

MAR11280 

* 

MAR11290 

IF(LOS  .LT.  SL0S1  .OR.  LOS  .  GT.  SL0S2)  GO  TO  1 

MAR11300 

IL=L0S-3L0S1+1 

MAR11310 

** 

YEARS  OVER  ST2  ARE  USED  FOR  VALIDATION 

MAR11320 

IF( YR  .LT.  ST1)  GO  TO  1 

MAR11330 

IT=YR-ST1+1 

MAR11340 

** 

COMPUTE  TOTAL  LOSS 

MAR11350 

Y=Q 

MAR11360 

DO  50  1=1,8 

MAR11370 

Y=Y  +  LOSS(I) 

MAR11380 

50 

CONTINUE 

MAR11390 

* 

MARI 1400 

RETURN 

MAR11410 

999 

IE0F=1 

MARI 1420 

END 

MAR11430 
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APPENDIX  C.  SIMULATION  FOR  CHOICE  OF  AVERAGE 
INVENTORY  VALUES 


SIMULATION  COMPARISON  OF  INVENTORY  MEAN  PERFORMANCE 

PURPOSE:  TO  IDENTIFY  THE  BEST  MEAN  VALUE  (ARITHMETIC, 

GEOMETRIC  AND  HARMONIC)  FOR  USE  IN  THE  INVERSION  OF 
TRANSFORMED  DATA  TO  THE  ORIGINAL  SPACE. 

PROGRAMMER:  CAPT  C  R  DICKINSON,  USMC 
VARIABLE  DESCRIPTION: 


PROGRAM  MCSIM 

iriHrirkirirkirirkfrMrkirMrk^iHrtrMrklrMrkllrlrklr*  ***  AA***'A***A»  A  A  iHHHhHWnHnlnHr*  A  AW**  A 
* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


LAMBDA 

TRIALS 

REPS 

BASIC 

FTE1 


PARAMETER  USED  IN  THE  POISSON  RANDOM  NUMBER 
GENERATOR 

VECTOR  OF  INVENTORY  VALUES  GENERATED  FROM 
THL  POISSON  RANDOM  NUMBER  GENERATOR 
NUMBER  OF  INVENTORY  VALUES  TO  BE  GENERATED 
BASIC  INVERSION  OF  TRANSFORMED  VALUES 
FREEMAN-TUKEY  EXACT  INVERSION  WITH 


FTE2  -  FREEMAN-TUKEY  EXACT  INVERSION  WITH 


Z 

MEANS 


★ 

* 

* 

★ 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


-  TRANSFORMED  SCALE  VARIABLE 

-  VECTOR  OF  ARITHMETIC,  GEOMETRIC  AND  HARMONIC  * 

* 
* 
* 
* 
* 
* 
* 
* 
* 
* 


MEANS 

TRSUM  -  SUMMATION  OF  TRIALS  FOR  USE  IN  THE 
ARITHMETIC  MEAN 

TRMULT  -  MULTIPLICATION  OF  TRIALS  FOR  USE  IN  THE 
GEOMETRIC  MEAN 

TRINV  -  SUMMATION  OF  THE  INVERSE  OF  TRIALS  FOR  USE 
IN  THE  HARMONIC  MEAN 

INPUT/OUTPUT: 

INPUT  -  HARD  CODED  PARAMETERS 
OUTPUT  -  SIM  OUT  A1  (DISK  FILE  03) 


*irMrirk*1i-/r-k1rk-frk-/i:MrirMr)rtrir/rir/r/firk-lrlrifA  iHnMW  *****  ******  *********  AnHUrftfrA1***** 


VARIABLE  DECLARATION 

INTEGER  TRIALS ( 100) ,Y( 100) ,REPS(3) ,PR0B,R, 

C  COUNT, TRSUM, YTEMP.YSUM 

REAL  Z( 100) ,T( 100) ,MEAN(3) ,BASIC( 3) ,DIFF(3) ,LAMBDA(5) , 
C  P(4),FTE1(3),FTE2(3) .TRMULT 

DOUBLE  PRECISION  DPSEED 

CALL  EXCMS( 'FILEDEF  02  DISK  REPLIC  LISTING  Al’) 


VARIABLE  INITIALIZATION 
PI=3. 14159 
DPSEED=889246.  DO 
LAM=5 
R=3 

PR0B=4 

C0UNT*2 


MCS00010 

MCS00020 

MCS00030 

MCS00040 

MCS00050 

MCS00060 

MCS00070 

MCS00080 

MCS00090 

MCS00100 

MCS00110 

MCS00120 

MCS00130 

MCS00140 

MCS00150 

MCS00160 

MCS00170 

MCS00180 

MCS00190 

MCS00200 

MCS00210 

MCS00220 

MCS00230 

MCS00240 

MCS00250 

MCS00260 

MCS00270 

MCS00280 

MCS00290 

MCS00300 

MCS00310 

MCS00320 

MCS00330 

MCS00340 

MCS00350 

MCS00360 

MCS00370 

MCS00380 

MCS00390 

MCS00400 

MCS00410 

MCS00420 

MCS00430 

MCS00440 

MCS00450 

MCS00460 

MCS00470 

MCS00480 

MCS00490 
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on  no  on 


* 

*** 


20 

* 


*** 


DATA  REPS/ 20, 40, 60/ 

DATA  LAMBDA/1.  ,2.  ,4. ,8.  ,16./ 
DATA  P/.  05,.  1,.  2,-4/ 


HEADER  PRINT 

WRITE  ( 02,’ C/T6,”  LAMBDA'  '  .T15,  "REPS'  '  ,T27,"P' 
T34,  "METHOD'  '  ,T46.'  'ARITHMETIC'  '  ,T59, 

' 'GEOMETRIC' ' ,T71, ' ’HARMONIC" /)' ) 


LOOP  FOR  LAMBDA  PARAMETER 
500  1=1, LAM 
DO  450  J=1,R 

CALL  GGP0S( LAMBDA( I ) , DPSEED , REPS( J) , TRIALS , IER) 
NSUM=0 

DO  10  J1=1,REPS( J) 

TRIALSC J1)=TRIALS( Jl)+1 
NSUM=NSUM+TRI ALS ( J 1 ) 


LOOP  FOR  PROBABILITY  PARAMETER 
DO  400  K=l,PROB 
C0UNT=C0UNT+4 

LINE  COUNTER  FOR  PAGE  BREAK 
IF  (COUNT.  GT.  55)  THEN 
WRITE  (2, ' ( 1H1) ' ) 

COUNT=2 

WRITE  (02, '  (/T6,  "LAMBDA’  '  ,T15,  ’  'REPS'  ’  ,T27, '  'P"  , 
T34,  "METHOD* '  ,T46  '  'ARITHMETIC'  '  ,T59, 
"GEOMETRIC'  ’  ,T71 , f ' HARMONIC" /)'  ) 

END  IF 


TRSUM=0 

TRINV=0 

TRMULT=0. 

ARCSUM=0. 

TSUM=0. 

YSUM=0 


LOOP  FOR  REPITITION  PARAMETER 
DO  20  Kl=l ,R£PS( J) 

CALL  GGBN( DPSEED, 1, TRIALS (Kl) ,P(K) ,YTEMP) 

YSUM=YSUM+YTEMP 

TRSUM=TRSUM+TR I ALS ( K 1 ) 

TRMULT=TRMULT+ALOG( REAL( TRI ALS( Kl ) ) ) 
TRINV=TRINV+1.  /REAL( TRIALS ( Kl) ) 

Y(K1)=YTEMP 

Z(K1)=.  5*(TRIALS(K1)+.  5)**.  5*(ASIN(2. *Y(K1)/ 
(TRIALS(K1)+1.  )-l.  )+ASIN(2.*(Y(Kl)+l.  )/ 
(TRIALS(K1)+1.  )-l. )+PI) 
T(K1)=Z(K1)/((TRIALS(K1)+.  5)**. 5) 
ARCSUM=ARCSUM+Z( Kl) 

TSUM=TSUM+T( Kl ) 

ZAVG=ARCSUM/REAL( REPS( J) ) 

TAVG=TSUM/REAL(REPS( J) ) 

MEAN  VALUE  CALCULATION 


MCS00500 
MCS005 10 
MCS00520 
MCS0Q53Q 
* MCS00540 
MCS00550 
MCS00560 
MCS00570 
MCS00580 
MCS00590 
MCS00600 
MCS00610 
MCS00620 
MCS00630 
MCS00640 
MCS00650 
MCS00660 
MCS00670 
MCS00680 
MCS00690 
MCS00700 
MCS00710 
MCS00720 
MCS00730 
MCS00740 
MCS00750 
MCS00760 
MCS00770 
MCS00780 
MCS00790 
MCS00800 
MCS00810 
MCS00820 
MCS00830 
MCS00840 
MCS00850 
MCS00860 
MCS00870 
MCS00880 
MCS00890 
MCS00900 
MCS00910 
MCS00920 
MCS00930 
MCS00940 
MCS00950 
MCS00960 
MCS00970 
MCS00980 
MCS00990 
MCS01000 
MCS01010 
MCS01020 
MCS01030 
MCS01040 
MCS01050 
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MEAN( 1 )=TRSUM/R£AL( REPS ( J) ) 

MEAN(2)=EXP(TRMULT*( 1.  /REAL(REPS( J)))) 

MEAN( 3 )=R£AL( REPS( J) ) /TRINV 

IF  (COS(TAVG).  GE.  0.  0)  THEN 
SIGNUM=1. 

ELSE 

SIGNUM=-1. 

ENDIF 

DO  30  K2=l ,3 

TAVG 1=ZAVG/ ( MEAN ( K2  ) +.  5)**. 5 
IF  (COS(TAVGl).GE. 0.0)  THEN 
SIG1-1. 

ELSE 

SIG1=-1. 

ENDIF 

BASIC  AND  FTE  INVERSE  METHOD  CALCULATION 

FTE1(K2)=. 5*( 1. -(SIGNUM*(1. -( ABS( SIN(TAVG)+ 

3  (SIN(TAVG)-(l.  /SIN(TAV3)))/MEAN(K2)))**2. )**. 5)) 

FTE2(K2)=.  5*( 1.  -(SIG1*(1. -( ABS( SIN(TAVG1)+ 

3  (SIN(TAVGl)-(l. /SIN(TAVG1)))/MEAN(K2)))**2. )**. 5)) 

BASIC(K2)=.  5*( 1.  +SIN( ZAVG/( (MEAN(K2)+.  5)**.  5) 

3  -(PI/2.))) 

OUTPUT  STATEMENTS 
DO  40  K3=l,3 

IF  (K3.EQ.  1)  THEN 

WRITE  (02, 1 (T6,F4.  1 ,T1S ,13 ,T25 ,F5. 3,T35 ,  " BASIC’ ' , 

:  T46 ,F7.  4,T59 ,F7. 4,T71 ,F7. 4) ' ) 

3  LAMBDA( I) ,REPS( J) ,P(K) , (BASIC( IK) , IK-1,3) 

ELSE 

IF  (K3.  EQ.  2)  THEN 

WRITE  (02, ' (T35 . ' 'FTEl' ' ,T46,F7.  4.T59.F7.  4, 

:  T71 ,F7.  4) ' )  (FTE1( IK) , IK=1 ,3) 

ELSE 

WRITE  (02,'(T35,"FTE2'  ’  ,T46,F7.  4,T59,F7.  4, 

3  T71.F7.4/)')  (FTE2(IK),IK=*1,3) 

ENDIF 

ENDIF 

CONTINUE 

CONTINUE 

CONTINUE 

CONTINUE 

STOP 

END 


MCS01060 
MCS01070 
MCS01080 
MCS01090 
MCS01100 
MCS01110 
MCS01120 
MCS01130 
MCS01140 
MCS01150 
MCS01160 
MCS01170 
MCS01180 
MCS0II9O 
MCS01200 
MCS01210 
MCS01220 
MCS01230 
MCS01240 
MCS01250 
MCS01260 
MCS01270 
MCS01280 
MCS01290 
MCS01300 
MCS01310 
MCS01320 
MCS01330 
MCS01340 
MCS01350 
MCS01360 
MCS01370 
MCS01380 
MCS01390 
MCS01400 
MCS01410 
MCS01420 
MCS01430 
MCS01440 
MCS01450 
MCS01460 
MCS01470 
MCS01480 
MCS01490 
MCS01500 
MC301510 
MCS01520 
MCS015 JO 
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